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PREFACE 

The  work  reported  herein  was  conducted  by  the  Arnold  Engineering  Development 
Center  (AEDC),  Air  Force  Systems  Command  (AFSC),  at  the  request  of  the  AEDC/DOT. 
The  AEDC  Project  Manager  was  Mr.  A.  F.  Money.  The  results  of  the  research  were 
obtained  by  ARO,  IncM  AEDC  Division  (a  Sverdrup  Corporation  Company),  operating 
contractor  for  the  AEDC,  AFSC,  Arnold  Air  Force  Station,  Tennessee,  under  ARO  Project 
Number  P32E-39D.  The  manuscript  was  submitted  for  publication  on  December  17,  1979. 

This  report  is  the  second  in  a  series  of  four  volumes  entitled  “Store  Separation  Testing 
Techniques  at  the  Arnold  Engineering Development  Center.”  Subtitles  of  these  volumes  are 
as  follows: 


Volume  I 
Volume  II 

Volume  III 

Volume  IV 


An  Overview 

Description  of  Captive  Trajectory  Store  Separation  Testing  in  the 
Aerodynamic  Wind  Tunnel  (4T) 

Description  and  Validation  of  Captive  Trajectory  Store  Separation 
Testing  in  the  von  Karman  Facility 

Description  of  Dynamic  Drop  Store  Separation  Testing 
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1.0  INTRODUCTION 

In  the  Arnold  Engineering  Development  Center  (AEDC)  Aerodynamic  Wind  Tunnel 
(4T),  store  separation  testing  with  the  Captive  Trajectory  Support  (CTS)  mechanism  was 
initiated  in  1968.  A  description  of  the  initial  CTS  hardware  and  the  separation  trajectory 
applications  program  are  contained  in  Ref.  1.  In  the  ensuing  years,  numerous  improvements 
to  the  operating  system  and  applications  program  have  been  implemented  as  a  result  of 
experience  accumulated  with  repeated  use  of  the  system.  The  purpose  of  this  report  is  to 
document  the  current  CTS  hardware,  the  systems  operation,  and  the  present  trajectory 
applications  program.  In  so  doing,  this  volume  will  automatically  supersede  Ref.  1.  This 
document  is  the  second  in  a  series  of  four  volumes  which  describe  store  separation 
capabilities  at  the  AEDC.  Volume  1  gives  an  overview  of  the  various  store  separation 
techniques,  Volume  III  describes  store  separation  testing  in  the  AEDC  Supersonic  Wind 
Tunnel  (A),  and  Volume  IV  covers  dynamic  drop  testing  capabilities  for  all  AEDC  wind 
tunnels. 


2.0  APPARATUS 


2.1  TEST  FACILITY 

The  Aerodynamic  Wind  Tunnel  (4T)  is  a  closed-loop  continuous  flow,  variable  density 
tunnel  in  which  the  Mach  number  can  be  varied  from  0.1  to  1.3  and  can  be  set  at  discrete 
Mach  numbers  of  1.6  and  2.0  by  placing  nozzle  inserts  over  the  permanent  sonic  nozzle.  At 
all  Mach  numbers,  the  stagnation  pressure  can  be  varied  from  300  to  3,700  psfa.  The  test 
section  is  4  ft  square  and  12.5  ft  long  with  perforated,  variable  porosity  (0.5-  to  10-percent 
open)  walls.  It  is  completely  enclosed  in  a  plenum  chamber  from  which  the  air  can  be 
evacuated,  allowing  part  of  the  tunnel  airflow  to  be  removed  through  the  perforated  walls  of 
the  test  section.  A  more  complete  description  of  the  test  facility  may  be  found  in  Ref.  2. 

During  captive  trajectory  testing,  two  separate  and  independent  support  systems  are 
used.  The  Captive  Trajectory  Support  (CTS)  and  the  wind  tunnel  main  pitch  sector  are  used 
to  support  the  store  and  aircraft  models,  respectively.  The  aircraft  model  is  supported  on  an 
adapter  sting  assembly  mounted  to  the  boom  of  the  main  pitch  sector. 

The  store  model  is  supported  on  a  sting  assembly  mounted  to  the  CTS  rig.  An  isometric 
drawing  of  a  typical  store  separation  installation  is  shown  in  Fig.  1,  along  with  a  block 
diagram  of  the  computer  control  loop  used  with  the  CTS.  A  schematic  showing  the  test 
section  details  and  the  location  of  typical  models  in  the  tunnel  is  shown  in  Fig.  2.  A 
photograph  showing  a  CTS  test  installation  with  multiple  positions  of  the  released  store  is 
shown  in  Fig.  3.  Further  description  of  the  CTS  rig  can  be  found  in  Ref.  2. 
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2.2  CAPTIVE  TRAJECTORY  SUPPORT  SYSTEM 

2.2.1  General 

The  CTS  is  used  primarily  for  the  trajectory  analysis  of  air-launched  stores  as  a 
separation  simulator  which  uses  the  wind  tunnel  as  a  six-degree-of-freedom  function 
generator  for  the  aerodynamic  coefficients  of  the  store.  The  CTS  hardware  consists  of  a  six- 
degree-of-freedom  store  model  support  with  a  closed-loop,  analog-control  positioning 
system  for  each  degree  of  freedom  and  interface  equipment  to  provide  communications  with 
the  AEDC  Propulsion  Wind  Tunnel  Facility  (PWT)  computer.  The  CTS  model  support  and 
positioning  systems  were  designed  and  built  by  General  Dynamics,  Convair  Division.  The 
interface  hardware  and  the  software  required  for  trajectory  generation  and  data  reduction 
were  developed  by  the  AEDC/PWT  Instrumentation  Branch. 

The  speed  and  precision  of  the  CTS  position  control  promote  its  use  for  nontrajectory 
tests  also.  A  sequence  of  positions  can  be  rapidly  traversed  with  the  desired  data  collected  at 
each  point.  Typical  uses  to  date  are  listed  below: 

1.  Grid  test:  the  store  model  is  located  at  various  positions  and  attitudes  relative  to 
the  aircraft  model;  forces  and  moments  are  measured,'  and  aerodynamic 
coefficients  are  calculated  and  displayed. 

2.  Free-air  test:  the  store  model,  with  no  aircraft  model  present,  is  rotated  in  pitch, 
yaw,  and/or  roll,  and  data  are  reduced  as  in  the  grid  test. 

3.  Flow-field  survey;  a  pressure  probe  mounted  on  the  CTS  rig  is  used  to  map  any 
region  of  interest. 

2.2.2  Slore  Model  Support 

The  CTS  is  an  electromechanical  system  with  six  degrees  of  freedom.  All  axes  of  motion 
are  contained  within  a  single  mechanism  that  is  independent  of  the  aircraft  model  support. 
Drive  motors  located  in  a  housing  attached  to  the  tunnel  structure  above  the  tunnel  diffuser 
are  printed  circuit  armature,  d-c  electric  motors  with  extremely  fast  response.  The  motors 
will  come  up  to  speed  in  approximately  O.J  sec.  The  motors  for  axial  and  vertical  motion  are 
rated  at  780  in.-oz  of  torque  at  a  speed  of  1,060  rpm.  For  pitch,  yaw,  roll,  and  transverse 
horizontal  motion,  the  motors  are  rated  at  120  in.-oz  of  torque  at  a  maximum  speed  of  2,750 
rpm.  The  horizontal,  pitch,  and  yaw  maximum  velocities  have  been  reduced  by  a  factor  of 
ten  to  minimize  overtravel  in  the  event  of  aircraft-store  fouling  (see  Section  2.2.5).  The 
resulting  linear  and  angular  velocities  of  the  six  degrees  of  freedom  are  as  follows: 
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Component 

Xr,  Axial 
Yr,  Horizontal 
Zr,  Vertical 

vr,  Pitch 
*?Ri  Yaw 
Ur,  Roll 


Velocity 

1.7  in. /sec 
0.5  in. /sec 

1.1  in. /sec  down, 

2.2  in. /sec  up 
2.0  deg/sec 
2.0  deg/ sec 

55.0  deg/ sec 


The  axial,  vertical,  and  horizontal  motions  are  accomplished  by  driving  bail  screws.  The 
envelope  of  translation  of  the  support  head  is  ±  15  in.  away  from  the  tunnel  centerline  in  the 
transverse  horizontal  and  vertical  directions.  The  axial  range  is  ±  18  in.  from  a  reference 
pitch  axis  location  at  tunnel  station  133.26.  Pitch  and  yaw  motions  are  accomplished  by 
driving  the  respective  gear  sector  with  a  conical  worm  gear  located  in  the  head  of  the 
support.  The  maximum  angular  range  of  pitch  and  yaw  motion  is  ±  45  deg.  Roll  motions  are 
accomplished  using  a  roll  shaft  driven  by  an  eccentric  gear  reduction  drive  with  a  maximum 
angular  motion  of  ±360  deg.  Zero,  3-in.,  and  6-in.  offset  roll  mechanisms  (Fig.  4)  are 
available  for  test  applications.  The  sting  support  and  balance  hardware  for  nonrolling  stings 
is  shown  in  Fig.  5. 

The  axial  and  vertical  motors  are  connected  to  their  respective  ball  screws  by  timing 
belts.  Power  for  transverse  horizontal  motion,  pitch,  yaw,  and  roll  is  transmitted  to  the 
support  head  by  flexible  shafts.  Position  readout  for  each  degree  of  freedom  is  accomplished 
by  the  use  of  precision  rotary  potentiometers  which  are  driven  with  a  minimum  gear 
reduction  between  the  motion  gear  and  the  potentiometer. 


2.2.3  Position  Control 


A  schematic  of  the  CTS  control  system  is  shown  in  Fig.  1.  Signal-conditioning  equipment 
and  position  control  and  monitoring  equipment  are  located  on  the  CTS  control  console 
shown  in  Fig.  6.  The  control  console  is  located  in  an  instrument  room  adjacent  to  the 
tunnel.  The  position  indicators  and  some  of  the  components  of  the  control  panel,  except  the 
manual  positioning  potentiometers  and  override  switches  visible  in  Fig.  6,  are  duplicated  on 
a  panel  of  the  Tunnel  4T  Data  Production  Console  (DPC)  for  monitoring  purposes.  The 
DPC  panel  also  contains  command  switches,  Fig.  7,  to  the  computer  for  initiating  and 
stopping  a  trajectory  and  controlling  output  trajectory  data  from  the  computer.  In 
computer-controlled  operation,  the  CTS  position  command  signals  are  applied  to  the 
summing  junctions  of  operational  amplifiers,  Fig.  1,  by  digital-to-analog  converters  (DAC) 
which  are  controlled  and  updated  by  the  computer.  For  manual  operation,  the  DAC  inputs 
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are  replaced  by  potentiometers  for  manually  positioning  the  CTS.  The  controllers  respond 
to  the  difference  between  commanded  and  actual  rig  positions  as  computed  by  the 
operational  amplifiers  and  drive  printed-circuit  motors  through  silicon-controlled  rectifier 
(SCR)  bridge  circuits.  Back-emf  of  the  motors,  sampled  when  the  SCR  bridges  are  turned 
off,  provides  velocity  feedback.  Motor  velocity  is  proportional  to  the  position  error  with 
maximum  speed  obtained  for  a  4-percent  error.  The  threshold  for  movement  corresponds  to 
an  error  signal  of  less  than  0.05  percent.  The  controllers  provide  motor  overload  protection 
by  electronically  limiting  the  drive  currents.  For  any  axis,  the  overall  positioning  error 
including  effects  of  rig  misalignment,  potentiometer  nonlinearity,  backlash,  power  supply 
drift,  and  other  error  sources  is  less  than  0.2  percent  of  full-scale  travel. 

2.2.4  Compuler-CTS  System  Interface 

A  general  block  diagram  of  the  CTS  system  and  computer  interface  is  shown  in  Fig.  1.  A 
strain-gage  balance  located  inside  the  store  model  detects  the  aerodynamic  forces  and 
moments  on  the  model.  The  resulting  force  and  moment  signals  are  processed  by  the  Digital 
Data  Acquisition  System  (DDAS).  The  Digital  Multiplexer  and  Control  System  computer 
(DMACS)  obtains  the  tunnel  conditions  and  the  six  CTS  positions.  The  CTS  positions, 
aerodynamic  data,  and  tunnel  conditions  are  then  input  to  the  facility  computer  which 
performs  the  prescribed  trajectory  calculations  and  concludes  the  cycle  by  transmitting  the 
six  new  calculated  positions  to  the  six  respective  digitai-to-analog  converters  (DAC).  At  the 
conclusion  of  rig  movement,  another  data  cycle  is  automatically  initiated. 

2.2.5  System  Safety  Provisions 

For  overtravel  protection,  a  dual  limit  switch  is  provided  in  each  direction  for  each 
degree  of  freedom.  If  the  first  limit  is  exceeded,  a  controller  safety  circuit  is  activated  which 
applies  dynamic  braking  to  stop  the  motor.  This  limit  can  be  overridden  to  return  to  a  safe 
position.  If  a  backup  limit  is  exceeded,  the  motor  controller  power  is  shut  off,  and  the  motor 
must  be  manually  cranked  back  into  the  operating  range. 

The  CTS  system  is  also  electrically  connected  to  automatically  stop  the  CTS  movement  if 
the  store  model  or  CTS  contacts  the  aircraft  model,  the  aircraft  support  sting,  or  the  test 
section  walls.  Television  monitors  located  in  the  Tunnel  4T  instrument  room  and  on  the 
Tunnel  4T  control  console  are  used  for  visual  observation  during  movement  of  the  store  to 
the  starting  position  and  during  the  controlled  positioning  of  the  store  in  the  trajectory. 
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Additional  protection  is  provided  by  mechanical  brakes  which  are  actuated  to  prevent 
horizontal,  vertical,  or  axial  movement  in  case  of  a  power  failure.  Brakes  are  not  required 
for  the  other  degrees  of  freedom  since  the  loads  are  not  large  enough  to  overcome  the 
friction  of  the  drive  train. 

2.2.6  Store  Model  Alignment  System 

In  Tunnel  4T,  the  pylons  and  racks  of  the  aircraft  models  contain  an  optical  sensor 
which  enables  the  store  model  to  be  accurately  positioned  at  the  carriage  position.  The 
optical  sensor  emits  infrared  radiation  and  detects  the  reflected  radiation  from  the  store 
model.  The  signal  of  the  reflected  radiation  is  inversely  proportional  to  the  distance  between 
the  store  and  the  sensor.  The  sensor  is  sensitive  to  store  positions  both  vertically  and  laterally 
with  respect  to  the  pylon  surface.  In  test  peculiar  cases,  the  pylons  and  racks  may  be 
instrumented  with  spring-loaded  plungers  (touch  wires)  which  are  electrically  connected  to 
give  a  visual  indication  on  the  DPC  and  control  console  when  the  store  model  makes  contact 
with  the  touch  wire. 

2.2.1  Store  Mode) 'Carriage-Positioning  Modes 

There  are  four  basic  operational  modes  for  positioning  the  store  model  for  initiating  a 
trajectory.  Initially,  the  store  is  manually  positioned  offline  in  its  carriage  position  and  the 
six  coordinates  (plus  aircraft  angle  of  attack)  constituting  a  “touch  point”  are  recorded  by 
the  facility  computer.  The  CTS  operator  located  at  the  DPC  selects  one  of  the  four 
operational  modes  for  positioning  the  store  to  its  initial  position. 

The  first  mode  allows  the  store  to  be  positioned  from  a  safekeeping  location  (nominally  I 
in.  away,  vertically)  by  making  a  series  of  small  movements  (sting  deflection  corrections 
included)  as  the  store  model  is  automatically  driven  toward  the  calibrated  position  by  the 
facility  computer.  Movement  is  terminated  when  the  proper  signal  is  sensed  by  the  optical 
sensor  or  the  touch  wire.  The  second  mode  of  operation  is  the  same  as  the  first  except  that 
the  store  model  is  moved  to  the  calibrated  spatial  coordinates  rather  than  to  the  preselected 
signal  from  the  optical  sensor  or  touch  wire.  The  third  operational  mode  allows  the  store 
model  to  be  moved  to  its  initial  position  in  one  motion  (including  the  necessary  sting 
deflections).  The  fourth  and  final  operational  mode  allows  the  CTS  operator  to  manually 
move  the  store  model  to  the  desired  spatial  coordinates  or  carriage  position. 


II 


AEDC-TR-79-1 


3.0  TRAJECTORY  GENERATION  PROGRAM 


3.1  GENERAL 

The  Tunnel  4T  trajectory  generation  applications  program  can  be  divided  into  three 
basic  blocks:  the  open-loop  service  routine  job,  the  CTS  rig  control  job,  and  the  trajectory 
generation  job.  The  open-loop  service  routine  can  be  accessed  by  the  data  acquisition  panel 
(DAP)  operator  at  many  points  and  contains  all  the  setup  and  calibration  programs 
necessary  for  the  conduct  of  the  test.  These  include  recording  pre-  or  postrun  instrument 
readings,  loading  program  constants,  initiating  the  summary  program,  performing  various 
calibrations,  performing  CTS  and  balance  instrumentation  checks,  obtaining  store  model 
weight  tare  values,  recording  model  check  loading  values,  and  recording  touch  points. 

The  CTS  rig  control  job  is  accessed  by  the  DAP  operator  primarily  in  three  ways.  First, 
“Initialize  Control”  is  performed;  this  transfers  control  of  the  CTS  rig  from  manual  to 
computer-controlled  operation.  Once  the  CTS  is  in  computer  control,  the  “Reset  CTS”  or 
“Set  CTS  to  IP”  function  may  be  activated  by  the  DAP  operator.  When  the  “Reset  CTS” 
function  is  chosen,  the  CTS  rig  is  automatically  moved  from  its  present  location  to  the 
station-keeping  position,  and  the  desired  aircraft  model  angle  of  attack  is  set.  When  the  “Set 
CTS  to  IP”  function  is  used,  the  store  model  is  automatically  moved  from  the  station- 
keeping  position  to  the  initial  position  (normally  carriage),  accounting  for  sting  deflections 
under  load  on  both  store  model  position  and  attitude.  The  station-keeping  position,  aircraft 
model  angle  of  attack,  and  store  model  initial  position  are  all  described  by  the  selected  touch 
point.  After  the  store  initial  position  is  determined  to  be  within  allowable  tolerances,  the 
closed-loop  trajectory  generation  process  (see  flow  diagram,  Fig.  8)  is  automatically 
activated.  The  new  CTS  coordinates  and  attitudes,  which  are  generated  from  the  trajectory 
equations,  are  set  by  the  closed-loop  module  of  the  CTS  rig  control  job.  The  closed-loop 
operation  will  be  continued  until  an  internally  defined  limit  is  reached  or  until  external 
termination.  The  closed-loop  process  may  be  paused  and  then  continued  or  terminated  by 
the  DAP  operator  at  any  time.  When  a  store  ground,  balance  static  limit,  balance  dynamic 
limit,  emergency  stop,  or  other  limit  is  encountered,  an  error  flag  is  set  which  automatically 
initiates  the  return  to  open-loop  operation.  A  detailed  description  of  the  CTS  mechanism  is 
given  in  Section  2,  and  details  of  the  trajectory  generation  process  are  contained  in  Sections 

3.2  through  3.7  and  Appendixes  A  through  L. 

The  general  composition  of  the  trajectory  generation  package,  as  outlined  in  the  flow 
diagram  of  Fig.  8,  consists  of  staging/initialization,  input  processing,  full-scale  trajectory 
calculations,  integration/extrapolation,  output  processing,  and  closed-loop  CTS 
positioning.  When  wind  tunnel  aerodynamic  data  are  utilized  in  the  full-scale  trajectory 
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calculations,  the  store  model  measured  forces  and  moments  are  reduced  to  coefficient  form 
and  applied  with  the  proper  full-scale  store  dimensions  and  flight  dynamic  pressure.  The 
equations  of  motion  allow  for  six-degrees-of-freedom  movement  of  the  released  store  and 
are  given  in  Appendix  I  (see  Fig.  9  for  definition  of  the  body-axis  system).  In  addition  to  free 
motion  releases,  the  equations  include  provisions  for  several  modes  of  staged  separation, 
aircraft  accelerated  flight  (Fig.  10),  and  aircraft  dive  or  bank  maneuvers.  Assumptions  and 
techniques  used  in  the  development  of  the  motion  equations  are  described  in  Volume  1  of 
this  series.  The  full-scale  force  and  moment  equations  (Appendix  H)  include  terms  to 
account  for  aerodynamic  damping,  weight,  thrust,  ejectorforces,  static  aerodynamic  forces, 
and  external  input  forces.  Integration  of  the  accelerations  and  velocities  is  accomplished 
using  the  Adams-Moulton  (predictor-corrector)  algorithm  with  a  Runge-Kutta  algorithm  to 
start  the  process. 

The  modular  structure  of  the  program  allows  for  the  implementation  of  new  routines 
that  might  be  required  for  a  specific  test  (e.g.,  autopilot  simulation)  with  a  minimum  of 
effort.  Requirements  for  the  addition  of  such  routines  and  for  standard  program  utilization 
are  given  in  Section  5. 

3.2  STAGING/INITIALIZATION 

The  staging  process  of  the  trajectory  generation  package  is  accomplished  in  two  steps. 
The  first  step  is  executed  by  a  constants  point  request  during  open-loop  operation.  The 
program  assembles  an  array  of  360  constants  from  previously  stored  permanent  files; 
performs  calculations  such  as  model  reference  areas  and  lengths,  transfer  distances,  CTS  rig 
physical  parameters,  and  full-scale  initial  position  coordinates  and  attitudes;  tabulates  the 
constants  and  calculated  parameters;  and  stores  the  constants  and  calculated  parameters  in  a 
temporary  initialization  disk  file.  Details  of  the  assembly  process  and  the  calculation 
equations  are  contained  in  Appendix  A.  The  second  step  of  the.  staging  process  is  executed  in 
the  first-pass  loop  of  the  trajectory  generation  job  (see  Fig.  8)  and  consists  of  reading  the 
temporary  initialization  and  CTS  rig  position  files  and  storing  the  values  in  the  common 
array  of  the  closed-loop  program. 


Initialization  of  the  trajectory  generation  job  is  also  carried  out  in  the  first-pass  loop  of 
the  closed-loop  program.  It  includes  assigning  and  initializing  the  integrators,  initializing 
program  control  flags  and  counters,  and  doing  once-only  calculations.  The  standard 
trajectory  package  uses  only  21  of  the  50  available  integrators,  so  the  remainder  may  be  used 
in  test  peculiar  applications.  The  once-onty  calculation  equations  are  listed  in  Appendix  A. 
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3.3  INPUT  PROCESSING 

The  primary  functions  of  the  input  processing  job  are  to  acquire  the  experimental  data, 
define  wind  tunnel  test  conditions,  calculate  the  measured  aerodynamic  coefficient  data, 
determine  angular  and  linear  deflections  of  the  model  balance/sting  combination  under 
load,  and  manage  and  update  the  extrapolation  data  base.  Secondary  functions  include  data 
acquisition  and  tunnel  conditions  validity  checks,  deflection  change  limit  checks,  and  critical 
moment  checks  on  the  CTS  sting  and/or  gears.  Sequencing  of  the  input  processing  events  is 
described  in  the  flow  diagram  of  Fig.  8. 

The  wind  tunnel  test  conditions  are  calculated  using  a  standard  Tunnel  4T  data  reduction 
routine.  Balance  readings  are  converted  into  gross  forces  and  moments  using  a  standard 
PWT  six-component  balance  data  reduction  program.  The  equations  which  are  used  to 
calculate  the  aerodynamic  coefficients  and  sting  deflections  from  the  gross  forces  and 
moments  are  given  in  Appendix  B.  If  differences  in  sting  deflections  between  succeeding 
points  exceed  allowable  tolerances  (0.05  in.  for  linear  deflections,  0.15  deg  for  pitch  and 
yaw,  and  1.0  deg  for  roll),  subsequent  calculations  are  bypassed  and  the  model  is 
repositioned  using  the  updated  deflections.  This  check  was  incorporated  because  model 
positioning  is  always  based  on  the  sting  deflections  calculated  for  the  previous  loading 
condition. 

For  critical  moment  checks,  the  model  gross  forces  and  moments  are  resolved  into 
moment  loadings  about  the  CTS  pitch  and  yaw  gears.  If  these  moments  exceed  the 
maximum  allowable  of  900  in. -lb,  the  “critical  moment”  flag  is  set  and  the  trajectory  is 
terminated  by  the  output  processing  job  before  another  rig  movement  is  executed.  If  either 
the  zero  or  the  6-in,  offset  roll  mechanism  is  used  in  the  balance/sting  combination,  a  more 
stringent  moment  check  is  required  which  is  described  in  the  Test  Facilities  Handbook  (Ref. 
2).  CTS  sting  or  gear  moment  overloads  are  not  often  encountered  during  testing  since  the 
static  load  capacities  of  most  CTS  balances  are  reached  well  before  critical  moments  would 
be  obtained. 

On  each  pass  through  input  processing,  the  extrapolated  coefficient  Values  are  identified 
and  set  equal  to  the  current  measured  values  (see  Appendix  B)  so  that  the  final  evaluation  of 
the  trajectory  equations  (before  data  output)  is  based  on  measured  aerodynamic  data.  Two 
files  are  maintained  for  extrapolation  purposes,  and  each  file  contains  three  values  for  each 
of  the  identified  extrapolation  parameters.  The  tunnel  data  file  retains  the  current  and  two 
previously  measured  values  of  each  parameter  and  is  updated  only  during  input  processing 
such  that  (typically): 
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(Cn).  =  CN 

Therefore,  when  fully  established,  the  incremental  full-scale  time  spacing  of  the  tunnel  data 
file  corresponds  to  the  program  data  acquisition  time  increment  (i.e.,  values  are  stored  at 
t  =  0,  0.5At,  At,  1.5At,  2At,  3At,  4At,  .  .  .).  These  coefficient  stacks  are  used  to  generate 
pseudo  data  in  the  multiple-pass  integration  routine  (see  Section  3.5). 

The  extrapolator  file  is  updated  on  each  pass  through  input  processing  but  is  also 
updated  by  the  pseudo  data  generated  in  the  multiple-pass  integration.  When  fully 
established,  incremental  full-scale  time  spacing  of  the  extrapolator  file  corresponds  to  the 
program  integration  time  increment  (i.e.,  values  are  stored  at  t  -  0, 0.55t,  fit,  1.55t,  2fit,  35t, 
46t  .  .  .).  One  absolute  and  two  incremental  coefficient  values  are  stored  such  that 
(typically): 

(“»)„  -  <4Ck),-. 

(“«)„  -c...- (=■-.), 

(CN,x).  =  CN.X 

The  incremental  coefficient  stacks  are  used  for  extrapolation  in  the  corrector  loop  of  the 
integration  algorithm.  The  basic  program  uses  6  of  the  25  available  extrapolators;  the 
remainder  may  be  used  in  test  peculiar  applications. 

3.4  TRAJECTORY  CALCULATIONS 

3.4.1  General 

In  the  trajectory  calculations  routine,  the  sequence  of  operations  is  to  update  the 
integration  parameter  buffers,  evaluate  the  trajectory  generation  equations,  and  then  update 
the  integrator  buffers.  The  trajectory  generation  equations  are  divided  into  functional 
modules  to  facilitate  handling.  The  equations  and  flow  diagrams  for  each  of  the  functional 
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modules  are  contained  in  Appendixes  C  through  I.  A  brief  description  of  each  module  in  the 
sequence  of  evaluation  is  given  as  follows. 

3.4.2  Conversion  Module  (Appendix  C) 

Several  diverse  functions  are  performed  in  the  conversion  module,  including 
assigning  the  appropriate  positions  and  velocities  to  the  integrator  results,  defining  the 
inertial-to-body-axis  direction  cosine  matrix,  describing  the  store  position  and  linear  velocity 
with  respect  to  the  inertial-axis  system  origin,  and  describing  the  store  attitude  with  respect 
to  the  inertial-axis  system  coordinate  directions.  The  nonrolling  sting  restraint  is 
implemented  in  this  module  when  required  (see  Appendix  L  for  information  concerning  this 
restriction).  Additional  store-related  calculations  include  body-axis  weight  components, 
total  velocity  (with  respect  to  a  space-fixed  axis  system),  dynamic  pressure  at  altitude,,  angle 
of  attack,  sideslip  angle,  and  lanyard  length. 

3.4.3  Offset  Coefficient  Module  (Appendix  D) 

The  purpose  of  this  module  is  to  define  any  aerodynamic  data  inputs  which  are  required 
in  addition  to  the  measured  aerodynamic  coefficients.  The  standard  offset  coefficient 
module  allows  for  constant  coefficient  inputs  (as  might  typically  be  used  to  correct  for 
reduced-scale  model  asymmetries)  and  for  a  ramp  axial-force  input  (to  simulate  drogue 
chute  deployment).  On  a  test  peculiar  basis,  specialized  modules  can  be  routinely  substituted 
for  the  standard  module  when  more  complex  offset  aerodynamic  coefficient  inputs  are 
required  (e.g.,  autopilot  simulation). 

3.4.4  Total  Coefficient  Module  (Appendix  E) 

This  module  sums  measured,  offset,  and  aerodynamic  damping  coefficient 
contributions. 

3.4.5  Thrust  Module  (Appendix  F) 

The  thrust  module  defines  the  values  of  simulated  full-scale  thrust  forces  and  moments. 
For  the  standard  module,  one  or  two  fifth-degree  polynomial  curve  fits  are  used  to  describe 
the  longitudinal  thrust  force  as  a  function  of  time.  Only  longitudinal  thrust  forces  are 
considered,  and  moments  arise  only  through  jet-damping  coefficient  contributions.  Module 
options  include  thrust  force  simulation  with  no  delay,  onset  of  thrust  forces  delayed  by  a 
time  increment,  and  onset  of  thrust  forces  delayed  by  a  lanyard  pull  followed  by  a  time 
increment  (see  Figs.  F-l  and  F-2).  In  test  peculiar  applications,  thrust  vector  control  could 
be  included  in  this  module. 
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3.4.6  Ejector  Module  (Appendix  G) 

The  ejector  module  defines  the  values  of  simulated  full-scale  ejector  forces  and 
moments.  For  the  standard  module,  two  independent  ejector  forces  may  be  input  using  one 
or  two  fifth-degree  polynomial  curve  fits  for  each  to  describe  the  forces  as  functions  of  time 
or  displacement  (see  Fig.  G-2).  Moment  contributions  are  determined  from  the  relationships 
of  the  ejector  piston  locations  with  respect  to  the  store  center  of  gravity,  and  duration  of 
ejector  action  may  be  defined  in  terms  of  either  time  or  displacement.  Module  options 
include  the  following  combinations:  ejector  forces  and  cutoff  =  f(t);  ejector  forces  and 
cutoff  =  ^displacement);  or  ejector  forces  =  f(t)  with  cutoff  =  ^displacement). 

3-4.7  Full-Scale  Force  and  Moment  Module  (Appendix  H) 

The  full-scale  forces  and  moments  resulting  from  weight,  aerodynamic,  thrust,  and 
ejector  contributions  are  summed  in  this  module. 

3.4.8  Dynamic  Equations  of  Motion  Module  (Appendix  1) 

In  addition  to  unrestrained  motion,  this  module  allows  for  the  following  three  basic 
types  of  staged  separation:  pivot  motion,  rail  launch  motion,  and  ejector-plane  motion  (see 
Fig.  1-3).  For  pivot  and  rail  motion,  the  store  is  initially  constrained  to  rotate  about  a  point 
other  than  the  center  of  gravity.  For  the  pivot  case,  the  rotation  center  is  fixed  with  respect 
to  the  aircraft  (a  typical  fuel  tank  release),  while  in  the  rail  case,  the  rotation  center  (hook)  is 
allowed  to  translate  with  respect  to  the  aircraft  (a  typical  missile  rail  launch).  For  ejector- 
plane  motion,  the  store  rotates  about  the  center  of  gravity  but  is  constrained  to  translate  and 
rotate  only  in  the  plane  of  the  ejectors  for  the  duration  of  ejector  action. 

The  equations  of  motion  for  each  option  are  solved  to  determine  the  linear  and  angular 
body-axis  accelerations,  the  inertial-axis  linear  velocities,  the  inertial-to -body-axis  direction 
cosine  derivative  matrix,  and  the  hook  accelerations  for  staged  release.  Each  of  these  terms 
is  then  assigned  to  the  appropriate  integrator  input.  Additional  module  output  includes 
reaction  forces  and  moments  required  to  impose  the  staged  separation.  Other  restrained 
motion  applications  could  be  incorporated  for  test  peculiar  requirements, 

3.5  INTEGRATION  AND  EXTRAPOLATION 

Integration  of  the  store  velocities  and  accelerations  is  accomplished  using  an  Adams- 
Moulton  algorithm  with  a  Runge-Kutta  start  (Ref.  3  and  Fig.  1 1).  The  Runge-Kutta  process 
requires  four  evaluations  of  the  trajectory  equations  over  each  integration  step  and  is  used 
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for  the  first  three  integration  time  increments.  Since  the  first  four  passes  through  the 
integration  module  using  the  Runge-Kutta  calculations  advance  time  only  in  half  steps,  five 
passes  are  required  to  accomplish  the  first  three  integration  time  steps.  At  the  beginning  of 
the  fourth  integration  step,  the  required  derivative  history  files  have  been  established  and  the 
Adams-Moulton  procedure  is  initiated.  The  Adams-Moulton  process  is  a  predictor- 
corrector  method  which  requires  two  evaluations  of  the  trajectory  equations  for  each 
integration  time  step.  Details  of  both  numerical  integration  procedures  are  contained  in  the 
flow  diagram  (Fig.  11). 

To  permit  the  trajectory  calculations  in  the  corrector  loop  of  the  integration  process  to 
be  made  independently  of  data  acquisition,  extrapolated  values  of  the  aerodynamic 
coefficients  are  used  in  the  calculations  made  after  time  advance  (see  Fig.  11).  The 
coefficients  are  extrapolated  using  a  quadratic  fit  of  the  form  shown  in  Fig.  12.  The  different 
values  of  the  extrapolation  equation  constants  used  during  passes  4  and  5  result  from  the 
initial  uneven  time  spacing  of  the  coefficient  values  stored  in  the  extrapolator  file  (see 
Section  3,3). 

A  second  extrapolation  procedure  is  employed  in  the  program  to  allow  multiple  passes 
through  the  integration  module  for  each  data  acquisition  cycle  (see  flow  diagrams  of  Figs.  8 
and  13).  This  permits  data  acquisition  time  to  remain  large  (typically,  10  to  20  msec,  full- 
scale  time)  while  integration  time  is  small  (typically,  0.5  to  1  msec,  full-scale  time).  As  a 
result,  wind  tunnel  test  time  required  for  each  trajectory  and  errors  associated  with  step 
functions  in  the  integrated  parameters  (see  Section  3.8)  can  be  minimized.  The  small  penalty 
extracted  by  the  multiple-pass  integration  (approximately  20  msec  computational  time  per 
pass  through  the  integrators)  is  more  than  adequately  compensated  by  the  acquired 
advantages. 

The  extrapolator  for  multiple-pass  integration  operates  on  the  tunnel  data  file  using  a 
quadratic  fit  of  the  form  described  in  the  module  flow  diagram  (Fig.  13)  and  generates 
pseudo  data  which  are  used  to  update  the  extrapolator  file.  The  pseudo-data  generated  by 
this  routine  are  used  in  the  trajectory  calculations  and  the  integration  procedure  exactly  as  if 
they  had  been  measured.  The  extensive  checking  at  the  beginning  of  the  routine  is  done  to 
insure  compatibility  with  the  Runge-Kutta  integration  and  to  set  the  proper  time  structure  of 
the  tunnel  data  file.  Manipulation  of  extrapolation  equations  in  passes  4  and  5  results  from 
the  initial  uneven  time  spacing  of  the  tunnel  data  file  coefficient  values  (see  Section  3.3). 

3.6  OUTPUT  PROCESSING 

The  output  processing  job  is  executed  in  two  phases.  The  first  phase  is  initiated 
immediately  after  completion  of  the  trajectory  calculations  which  follow  input  processing;  it 
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consists  of  calculating  additional  trajectory  parameters,  storing  pertinent  trajectory 
information  in  a  768-word  engineering-unit-data  disk  file,  tabulating  hard  copy  data 
(optional),  displaying  selected  trajectory  parameters  on  an  alphanumeric  cathode  ray  tube 
screen,  and  increasing  the  data  cycle/integration  time  increments  as  required  (see  flow 
diagram,  Fig.  14). 

The  additional  calculations  are  performed  just  prior  to  data  output  to  increase  program 
efficiency  since  these  parameters  are  useful  in  interpreting  trajectory  results  but  are  not 
required  in  the  trajectory  generation  process.  Additional  information  calculated  includes  eg 
displacement  and  store  attitudes  relative  to  the  flight,  nonrotating  flight,  pylon,  nonrotating 
pylon,  aircraft,  and  earth  axis  system  coordinate  directions;  store  nose  and  tail 
displacement  and  store  attitudes  relative  to  the  flight-,  nonrotating  flight-,  pylon-, 
nonrotating  pylon-,  aircraft-,  and  earth-axis  system  coordinate  directions;  store  nose  and 
tail  displacements  parallel  to  the  flight-  and  pylon-axis  coordinate  directions;  hook 
displacement  parallel  to  the  pylon-axis  coordinate  directions;  and  aerodynamic  coefficients 
in  the  stability,  wind,  and  aeroballistic  axes.  The  standard  tabulated  summary  data  printout 
is  shown  in  Table  1,  but  any  of  the  information  stored  in  the  disk  file  can  be  tabulated  as 
required. 

The  flow  diagram  for  the  data  cycle/integration  interval  increase  module  is  shown  in  Fig. 
14.  Options  are  to  double  the  data  cycle  and  integration  time  increments  or  double  the  data 
cycle  time  increment  only  when  total  store  displacement  (lanyard  length)  reaches  a 
prescribed  value.  However,  no  increase  is  allowed  regardless  of  the  prescribed  displacement 
value  until  the  derivative,  extrapolator,  and  tunnel  data  history  files  have  been  fully 
established.  This  option  is  exercised  to  speed  up  the  trajectory  generation  process  after  the 
store  has  reached  a  position  in  the  flow  Held  where  aerodynamic  coefficient  gradients  are  not 
expected  to  be  large.  This  completes  the  first  phase  of  output  processing,  and  integration 
follows  (see  flow  diagram,  Fig.  8). 

After  integration  is  completed,  the  second  phase  of  output  processing  is  begun.  It 
includes  advancing  the  data  cycle  pass  counter,  checking  the  critical  moment  flag  (see 
Section  3.3)  and  terminating  the  trajectory  if  so  required,  and  calculating  the  store  model 
position  coordinates  and  attitudes  which  will  be  executed  by  the  CTS  rig.  These  coordinates 
are  store  model  eg  displacements  with  respect  to  the  origin  of  the  flight-axis  system  (see  Fig. 
10)  and  angular  displacements  relative  to  the  free-stream  wind  vector  for  a  pitch,  yaw,  roll 
movement  sequence  (including  corrections  for  induced  angles  resulting  from  vertical  and 
lateral  velocities  of  the  store  eg,  see  Fig.  15).  The  equations  are  given  in  Appendix  J. 
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3.7  CTS  CLOSED-LOOP  POSITIONING 

The  primary  function  of  the  CTS  closed-loop  positioning  module  is  to  locate  the  store 
model  at  the  new  set  of  trajectory  coordinates  and  attitudes  obtained  from  output 
processing.  Secondary  functions  include  travel  limit,  sidewall  clearance,  and  DAP  terminate 
request  checks.  The  equations  and  flow  diagrams  for  this  module  are  contained  in 
Appendix  K. 

Positioning  of  the  store  model  is  accomplished  using  an  absolute  tunnel  coordinate 
system,  the  origin  of  which  is  described  by  the  midpoint  of  travel  of  the  CTS  linear 
positions.  After  the  store  model  is  located  at  the  initial  position  (by  the  touch  job),  the 
carriage  coordinates  (XTpi0,  Ytp,0,  Ztp,0)  are  calculated  by  the  touch  job  and  stored.  The 
same  equations  as  described  in  Appendix  K  are  used  except  that  the  carriage  coordinates 
rather  than  rig  positions  are  the  unknown  quanities. 

On  the  first  pass  through  the  module,  the  carriage  coordinates  and  constant  box  inputs 
(AYC,  AZc)  are  read  from  the  touch  file.  These  values  remain  constant  for  the  duration  of  the 
trajectory.  The  new  set  of  rig  angular  positions  is  calculated  by  substracting  sting  deflection 
and  sting  bend  angle  contributions  from  the  absolute  store  attitude.  It  should  be  noted  that 
the  sting  deflection  angles  used  are  calculated  from  the.  previous  aerodynamic  load 
condition,  but  experience  has  shown  that  deflection  changes  from  point  to  point  are 
normally  small.  If  changes  should  become  large,  corrections  are  made  in  input  processing 
(see  Section  3.3).  The  new  set  of  rig  linear  coordinates  is  determined  according  to  the 
contributions  resulting  from  CTS  rig  geometry,  linear  sting  deflections,  the  new  trajectory 
coordinates,  and  constant  box  inputs.  The  constant  box  inputs  are  included  in  the 
calculations  as  a  means  for  correcting  air-off  carriage  coordinates  for  air-on  conditions 
when  a  touch  sensor  is  not  located  at  that  particular  launch  station. 


If  the  new  CTS  linear  positions  or 'angular  orientations  are  calculated  to  be  outside  the 
normal  operating  limits  listed  below,  the  trajectory  is  automatically  terminated  before  a  rig 
movement  is  executed. 


CTS  Position  Drives 

Design 

Travel  Limits 

Axial,  in. 

±18 

Horizontal,  in. 

±15 

Vertical,  in. 

±15 

Pitch,  deg 

±45 

Yaw,  deg 

±45 

Roll,  deg 

±360 

Normal  Operational 
Travel  Limits 

Position 

Positive 

Negative 

Tolerances 

16.9 

-17.7 

±0.05 

14.8 

-14.7 

±0.05 

14.7 

-14.8 

±0.05 

43.0 

-44.3 

±0.15 

44.5 

-44.1 

±0.15 

359.0 

-359.0 

±1.0 
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Similarly,  if  the  new  position  of  the  store  nose  is  calculated  to  be  within  2  in.  of  the  tunnel 
wall,  the  trajectory  is  aborted. 

After  the  rig  movement  has  been  executed,  differences  between  the  command  and  set 
positions  and  angles  are  calculated.  IF  these  differences  are  within  allowable  tolerances  (see 
preceding  table),  the  trajectory  generation  process  is  continued.  If  not,  the  rig  is 
repositioned  until  the  tolerance  criteria  are  met. 

3.8  VALIDATION 

Verification  of  the  trajectory  generation  applications  program  was  accomplished  in  four 
basic  steps.  The  first  consisted  of  thoroughly  checking  the  computer  code  to  insure  accuracy 
and  completeness.  During  the  second  step,  the  outputs  of  the  individual  program  modules 
were  verified,  when  feasible.  During  the  third  phase,  trajectory  simulations  were  performed 
independently  of  the  wind  tunnel  for  verification  of  the  combined  equations.  The  fourth 
step  consisted  of  incorporation  of  the  trajectory  generation  package  with  the  wind  tunnel 
and  CTS  hardware  and  software  systems  for  total  program  verification. 

Because  of  the  modular  structure  of  the  program,  a  considerable  portion  of  the 
verification  was  accomplished  during  the  second  step-  Of  particular  importance  was  the 
validation  of  the  integration  scheme.  Even  for  such  relatively  extreme  motions  as  rapidly 
divergent  (30-deg  motion  in  0.1  sec)  and  short  period  sinusoidal  (15-deg  amplitude  with 
0. 3-sec  period),  errors  in  the  trajectory  results  derived  from  the  integration  algorithm  were 
determined  to  be  at  least  an  order  of  magnitude  less  than  errors  expected  from  normal  input 
data  scatter.  As  would  be  expected,  the  integration  algorithm  cannot  precisely  respond  to 
step  functions  in  the  accelerations  (e.g.,  ejector  cutoff,  staged  separation  termination).  For 
the  worst  case,  step  functions  effectively  increased  (or  decreased)  the  integrated  parameter 
values  by  one-half  integration  time  increment  at  the  point  where  the  step  occurred. 
Therefore,  using  a  small  integration  time  increment  in  the  trajectory  calculations  (by 
multiple  pass  integration,  see  Section  3,5)  can  minimize  step  function  effects. 

After  the  preliminary  verification  was  completed,  numerous  analytical  trajectories  were 
calculated  to  check  out  the  combined  trajectory  equations.  Typical  trajectory  results  are 
given  in  Figs.  16  through  20.  For  simple  gravity  releases,  the  motion  of  the  store  was 
identical  to  exact  solutions  of  the  equations  of  motion  (Fig.  16).  For  the  ejector-augmented 
gravity  release  (Fig.  16a),  the  calculated  motion  differed  only  slightly  from  the  exact  solution 
because  the  small  integration  time  increment  used  (1  msec)  minimized  the  step  function 
effect. 

Although  exact  solutions  could  not  be  calculated  for  comparison  with  the  pivoting  and 
rail-staged  separation  trajectories  because  of  the  complexity  of  the  motion  equations, 
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quantitative  verification  was  established  (Figs.  17  and  18).  Trends  in  the  trajectory  data  were 
consistent  with  the  applied  restraints,  and  termination  of  the  staged  separations  occurred  as 
expected.  Identical  trajectories  could  be  obtained  with  equivalent  inputs,  and  mirror  image 
solutions  could  be  accomplished  with  sign  changes  on  appropriate  inputs.  The  ejector  plane 
restraint  was  demonstrated  to  be  valid  (Fig.  19),  and  comparisons  of  trajectory  results  for 
the  present  and  Ref.  1  programs  were  quite  good  (Fig.  20).  Similarly,  experimental 
trajectories  from  Tunnel  4T  and  Tunnel  A  compared  favorably  (Ref.  4). 

4.0  AERODYNAMIC  AND  FLOW-FIELD  GRID  TEST  APPLICATIONS 

Approximately  75  percent  of  the  CTS  grid  applications  program  is  common  to  the 
trajectory  applications  program.  Only  a  few  changes  are  necessary  in  the  open-loop  service 
routine,  and  none  are  required  in  the  rig  control  job  (see  Section  3.1).  In  the  closed-loop 
portion  of  the  program,  all  trajectory-related  calculations  are  deleted  and  simply  replaced  by 
the  grid-positioning  algorithm.  In  addition,  when  a  flow-field  probe  (Fig.  21)  is  substituted 
for  the  model/balance  combination,  the  aerodynamic  coefficient  calculations  are  replaced 
by  the  flow  angle  calculations. 

For  grid  tests,  the  store  model  or  flow-field  probe  is  positioned  in  the  aircraft  flow  field . 
(or  free  stream)  at  selected  locations  and  orientations  which  are  preprogrammed  into  the 
digital  computer.  The  grid  matrix  can  be  defined  relative  to  any  coordinate  system  (e.g., 
flight  axis,  pylon  axis,  aircraft  axis),  and  the  choice  of  origin  location  and  positive 
coordinate  directions  is  arbitrary.  Since  the  grid  program  was  developed  primarily  for  use  in 
aircraft  model-related  test  applications,  translational  position  parameters  are  normally 
output  in  full-scale  feet.  However,  both  model-scale  and  full-scale  positions  can  be  made 
available  in  other  dimensional  units. 

Aerodynamic  grid  coefficient  data  are  normally  calculated  in  the  body-axis  system,  but 
stability-axis,  wind-axis,  or  aeroballistic-axis  coefficients  can  be  made  available.  Inter¬ 
ference  coefficient  values  (flow-field  aerodynamics  minus  free-stream  aerodynamics)  may 
also  be  determined  by  means  of  an  offline  data  reduction  program.  Aerodynamic  flow 
angles  are  calculated  relative  to  the  probe  axis  and  the  free-stream  wind  vector.  Standard 
tabulated  summary  printouts  for  aerodynamic  and  flow-field  grid  data  are  given  in  Tables  2 
and  3,  respectively. 

5.0  REQUIREMENTS  FOR  TESTS  USING  THE  CTS 
5.1  MODEL  DESIGN  CONSIDERATIONS 

In  addition  to  normal  design  requirements,  two  special  considerations  are  necessary  in 
model  design  for  test  programs  which  utilize  the  captive  trajectory  support  mechanism. 
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First,  it  is  imperative  that  all  models  be  electrically  conductive  to  insure  proper  operation  of 
the  store  ground  system  (see  Section  2.2).  If  this  safety  system  were  to  be  defeated,  extensive 
damage  to  balance  or  test  hardware  could  result  in  the  event  of  a  store/aircraft  model 
collision.  Second,  sting-supported  store  models  should  be  designed  for  minimum  weight  and 
with  the  model  mass  center  located  near  the  balance  electrical  center.  This  requirement  is 
imposed  to  alleviate  model/balance  dynamic  loading  problems  encountered  during  some  test 
programs. 

5.2  TEST  CRITERIA 

Unless  specified  otherwise,  the  values  of  the  trajectory  input  parameters  compiled  in  this 
section  are  required  to  be  constant  throughout  a  trajectory.  However,  for  specific  test 
applications,  selected  parameters  could  be  reprogrammed  as  variables  if  the  functional 
relationships  were  defined.  Values  of  the  following  parameters  (as  required)  must  be 
defined  in  the  test  planning: 

Store  Physical  Parameters  (Trajectory  or  Aerodynamic  Grid) 

K  A,  f|,  (2.  fj,  f,  Xcgl  Ycg,  Zcg,  hook  locations 
Grid  Matrix  Information  (Aerodynamic  or  Flow-Field  Grid) 


Origin  Location 

Orientations  and  positive  directions  of  grid  coordinates 
Orientation  of  store  with  respect  to  grid  coordinates 
Definition  of  grid  points 

Store  Mass  Properties  (Trajectory) 


Wt,  Ixxi  Ixy.  Ixz.  Iyy.  Iyz-  Izz 
Store  Aerodynamic  Inputs  (Trajectory) 

Qpt  Cmq,  C„r,  AXmiCgf  AXn<cg,  additional  aerodynamic  coefficient  inputs  (Ca.o),  etc.) 
Trajectory  Simulation  Parameters 


h,  Nz,  7,  <f>A/c 
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Ejector  Simulation  Parameters  (Trajectory) 

XFE,  AXaE.  Wm»  Z.HI,  Ze2 
(FE],  FE2)  versus  time  (stroke) 

Thrust  Simulation  Parameters  (Trajectory) 


tD»  ZE,  Fj.x  versus  t,  Cjdf,  Cjdm,  Cjdn 
Staged  Separation  Parameters  (Trajectory) 

X0,  Y0,  Zo,  Xpj,  Xpt2,  A0r,  staged  separation  mode 
Initial  Conditions  (Trajectory) 

For  trajectories  which  are  initiated  at  points  other  than  carriage,  the  initial  trajectory 
time  and  store  initial  positions,  orientations,  linear  velocities,  and  angular  velocities  must  be 
defined. 

Miscellaneous  Information  (Trajectory) 

If  lanyard  length  calculations  are  required,  attachment  coordinates  of  the  lanyard  to 
both  store  and  aircraft  should  be  defined. 

Autopilot  Applications  (Trajectory) 

The  simulation  of  trajectories  with  active  guidance  and  control  systems  requires  a 
mathematical  model  of  the  inertial  and/or  mechanical  response  of  the  systems.  Since  the 
mechanisms  are  unique  to  each  missile,  no  standard  programming  exists  to  describe  them. 
However,  the  standard  trajectory  program  is  capable  of  dealing  with  the  active  control 
situation  by  calculating  incremental  aerodynamic  coefficients  resulting  from  the  control 
surface  deflections.  Information  required  includes  a  mathematical  algorithm  describing  the 
control  surface  movements  as  functions  of  missile  acceleration,  velocity,  position,  attitude, 
etc.,  and  the  body-axis  aerodynamic  coefficient  variations  resulting  from  the  control  surface 
deflections.  Since  this  requires  test-unique  program  additions,  at  least  eight  weeks’  lead  time 
should  be  allowed  to  permit  program  preparation  and  checkout.  Sample  check  calculations 
should  be  provided,  if  available. 
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5.3  NOMENCLATURE  UPDATE 

Comparison  of  nomenclature  from  Volume  III  of  this  series  and  the  present  report  will 
disclose  differences  in  the  engineering  symbols  for  several  terms.  Since  the  publication  of 
Volume  111,  an  extensive  effort  has  been  initiated  at  AEDC  to  standardize  nomenclature 
among  the  test  facilities  and  among  different  types  of  test  programs,  and  differences  are 
generally  a  result  of  this  effort.  Where  terminology  differences  are  noted  between  the  two 
reports,  those  given  in  this  report  will  supersede  Volume  III  as  the  correct  notation. 

REFERENCES 

1.  Christopher,  J.  P.  and  Carieton,  W.  E.  “Captive-Trajectory  Store-Separation  System  of 
the  AEDC-PWT  4-Foot  Transonic  Tunnel.”  AEDC-TR-68-200  (AD839743),  September 
1968. 

2.  Test  Facilities  Handbook  (Eleventh  Edition).  “Propulsion  Wind  Tunnel  Facility,. 
Vol.  4.”  Arnold  Engineering  Development  Center,  June  1979. 

3.  Henrici,  Peter.  Discrete  Variable  Methods  in  Ordinary  Differential  Equations .  John 
Wiley  and  Sons,  Inc.,  New  York,  1964. 

4.  Hill,  D.  W.,  Jr.,  Best,  J.  T.,  and  Tolbert,  R.  H.  “Comparison  of  Store  Trajectory  and 
Aerodynamic  Loads,  and  Model  Flow-Field  Characteristics  Obtained  in  the  AEDC 
PWT/4T  and  VKF/A  Wind  Tunnels  at  Mach  Number  1.63.”  AEDC-TR -78-45 
(ADA065137),  February  1979. 


25 


Figure  1.  Isometric  drawing  of  a  typical  store  separation  installation 
and  a  block  diagram  of  the  computer  control  loop. 


TYPICAL  PERFORATED  WALL  CROSS  SECTION  TUNNEL  STATIONS  ANO 

DIMENSIONS  IN  INCHES 


STA.  STA.  STA.  STA. 

0.0  36.0  76.8  150.0 


Figure  2.  Schematic  of  the  tunnel  test  section  showing 
typical  model  locations! 
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Figure  3.  CTS  installation  photograph  showing  multiple 
locations  of  the  released  store. 
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DIMENSIONS  IN  INCHES 


b.  3-in.  offset 


BALANCE 


NAME 

A 

Q 16" 6.5  STRAIGHT  BALANCE 

5-.I6-.0065-.400M 

7.96 

0.16* 6.51b  SPEC  BALANCE 

5- .16 -.0065- SPEC 

nsEid 

O.IBB*  4.01b  SPEC  BALANCE 

9-  .186-. 0040-  SPEC 

1 1.76** 

0.30'  7.21b  SPEC  BALANCE 

6-. 30  -  .00  72- SPEC 

7.77** 

0.40"  10  lb  STRAIGHT  BALANCE 

6-. 40-  .010-. 40  M 

8.0  3 

0.40' 20>b  STRAIGHT  BALANCE 

6-.40-.020-.40M 

8.03 

0.40*  10  lb  SPEC  BALANCE 

6 -.40-  010-SPEC 

EXEhBi 

0.40"  201b  SPEC  BALANCE 

6-.40-. 020 -SPEC 

9.4  3** 

STING  SUPPORT 


L _ NAME _ 1 

B 

JTW-i  i;  M  MI  k  !  H  M  » 

ISi402F-8.00-.402M 

8.00 

1  1.25"  STRAIGHT  EXTENSION 

1 1.25 

13:50’  STRAIGHT  EXTENSION 

IS-.402 F -13.50- .402  M  | 

inm 

*GP  denotes  gage  point 
**  Includes  spec  sting 

d.  Additional  information 

Figure  4.  CTS  roll  mechanisms  and  associated  sting/balance  hardware. 
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YAW  PITCH 
CENTER  CENTER 


NAME 

NOMENCLATURE 

A 

0.16'  6.516  OFFSET  BALANCE 

4-.  16  -  .0065  -  3“  OFFSET-.  40  M 

1657 

0  40'  10  lb  OFFSET  BALANCE 

5 -  .40-  .010  -  3"  OFFSET- .40  M 

17.  62 

FOR  STING'  SUPPORT  INFORMATION  (B) 
SEE  FIGURE  4. 


a.  Offset  balance 


*GP  denote*  gage  point 


FOR  BALANCE  INFORMATION  (Al 
SEE  FIGURE  4. 


b.  Offset  stings 

Figure  5.  CTS  nonrolling  sting/balance  hardware. 
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.16"-  .30"  Bolonee  ond  CTS  Hardwore 

BALANCES 


NAME 

NOMENCLATURE 

■9 

.16*  6.5  lb  BALANCE 

5-  .IS -.0065- Spec 

6.64 

.IBB"  4.0  lb  BALANCE 

5-.  IBB -  0040-Spec 

664 

.30"  7.2  lb  BALANCE 

6-  .30-.0072.Spec 

2.65 

SPECIAL  STING  ADAPTER 


NAME 

NOMENCLATURE 

B 

2’  ADAPTER 

SPA  .30SSF-2.00-  30SSF 

2.00 

BENT  SPECIAL  STING 


NAME 

NOMENCLATURE  . 

12*  BENT  STING 

BSP.S. -Spec-6  635- 12*- .520  M 

B.58 

20*  BENT  STING 

BSP.S-Spee-B.T5l- 20*-. 520M 

8.52 

27*  BENT  STING 

BSP.S. -Spec-8.76  5 -27*-.520M 

8.43 

.40"  Balance  and  CTS  Hardware 


BALANCES 


NAME 

NOMENCLATURE 

A 

.40"  10  lb  BALANCE 

6-.40- .0(0-  Spec 

3.76 

.40'  20  lb  BALANCE 

6-. 40-  .020-  Spec 

3.76 

SPECIAL  STING  ADAPTER 


■DSH 

NOMENCLATURE 

B 

2"  ADAPTER 

SPA  .40SSP-2.00-.40  SSF 

2.00 

3.5"  ADAPTER 

SPA  .40 5SF-3.50  *  40  SSF 

3.50 

BENT  SPECIAL  ADAPTER 


NAME 

NOMENCLATURE 

C 

12*  BENT  STING 

BSP.S. -Spec  -9.30  -  12*-  .520 M 

6.9  5 

20*  BENT  STING 

BSP.S.-Sp*c -  9.23  -  20*-  .5 20 M 

867 

27*  BENT  STIM6 

BSP.S  -Spec -9. 14-27VS20M 

8.77 

BENT  STING  ADAPTER 


NAME 

NOMENCLATURE 

12*  ADAPTER 

BSA-.520F-  5. 1 25  -  12*9'  •  1.12  0  M 

4.38 

CTS  ADAPTER 


NAME 

NOMENCLATURE 

E 

3.5"  ADAPTER 

SA-I.I20F-  3.50-1. 1Z0M 

3.50 

10.0"  ADAPTER 

S  -  1.1  20F  -  10-00  -l.l 20  M 

IOlOO 

c.  Modular  stings  &  balances 
Figure  5.  Concluded. 
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:Null  Detectors 


Redundant 
Store  Ground 
Meter  Relay 


rrri 

% 

A 

a.  Control  console  components 
Figure  6.  Photographs  of  the  CTS  control  console. 
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Store  Ground 
and  Override 
Switch 


Position  Balance 
Indicator  Lights 


Redundant  Store 
Ground  Light 


Toxich  Wire 

Indicator 

Lights 


Manual  or 
Computer 
Control 
Switch 


Primary  and  Backup 
Limit  Indicators 


Manual  Positioning 
Potentiometers 


b.  Control  panel  components 
Figure  6.  Concluded. 
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Figure  8.  Flow  diagram  of  the  trajectory  generation 
closed-loop  program. 
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r  an 


FIRST  PASS 
(NPASS  =  0) 


INTEGRATION 
FLOW  DIAGRAM 


THIRD  PASS  SECOND/FOURTH  PASS  FIFTH  PASS 

(NPASS  *  2)  {NPASS-  I  OR  3)  (NPASS  •  4 ) 


SIXTH  PASS  — - 
(NPASS  2  5) 


OJ 

VO 


EXIT 


WHERE:  NPASS  IS  INTEGRATOR  PASS  COUNTER 

St  |jS  THE  INTEGRATION  TIME  INCREMENT 
P  IS  INTEGRATOR  OUTPUT 
P  IS  INTEGRATOR  INPUT 

PASS (1-5)  REPRESENTS 

RUNGE-KUTTA  INTEGRATION 
PASS (6""  )  REPRESENTS 

ADAMS-  MOULTON  INTEGRATION 

NOTE:  BETWEEN  EXIT  ONE  PASS 
AND  ENTER  NEXT  PASS, 
TRAJECTORY  CALCULATIONS 
ARE  PERFORMED. 


EXIT 


EXIT 


Figure  11.  Flow  diagram  for  the  integration  module. 


AEDC-TR-79 


AEDC-TR-79-1 


EXIT 

NOTE;  NXP  IS  THE  EXTRAPOLATOR  PASS  COUNTER 

Figure  12.  Flow  diagram  of  the  quadratic  extrapolation  routine. 
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Where  Kin,  .  NO  OF  I NTEGRAT  IONS / OATA  CYCLE  ( IF  *  I.  MUST  BE  EVEN  NO.) 
KPASS=  MULTIPLE  PASS  INTEGRATION  COUNTER 
PASS'  OATA  CYCLE  PASS  COUNTER 


Figure  13.  Flow  diagram  for  the  pseudo  data  generation  module. 
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STEP  OPTIONS 

0  NO  INCREASE 

1  DOUBLE  DATA  CYCLE  AND  INTEGRATION 

TIME  INCREMENT  FOR  ?LC  >  ESTEP 

2  DOUBLE  DATA  CYCLE  TIME  INCREMENT 

FOR  El£-  2 STEP 


EXIT 


NOTE!  NO  INCREASE  ALLOWED  UNTIL  THE  DERIVATIVE,  EXTRAPOLATOR,  AND 
TUNNEL  DATA  STACKS  HAVE  BEEN  FULLY  ESTABLISHED  AT  EVEN 
TIME  INCREMENTS 


Figure  14.  Flow  diagram  for  the  data  cycle/integration 
interval  increase  module. 
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OF  PITCH 


4Wi[»iPTi|liMliT« 


^1  1  I  I  I  I  I  I  I  I 


0  0.2  0.4  0.6  0.8  1.0 

t.sec 


b.  Angular  displacements 
Figure  16.  Concluded. 
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Figure  18.  Trajectory  results  for  rail-released  staged  separation, 
MOTION  =  7,  Wt  =  2,000  lb,  7  =  60  deg,  X0  =  5  ft. 


Table  1.  Standard  Trajectory  Tabulated  Summary  Data  Format 


RUN  TRAJ 

Mo. 

P. 
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00  Poo 

t_ 
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"^DP 

SH 

A 

h  8t 
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0,849  1200. 3  541.4  i/> 
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a 
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.0 
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1 
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0 
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X 

Y 

z 

e 
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Qs 
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Cm 
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cn 
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4a 
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Fe2 
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0.00 

0.  00 

0 . 98 

0 . 0 
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0.839 

0.212 
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3467 
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TO  .00 

-U‘,00 

0.02 

*U  .00 
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0,0 

1,22 

-0,00  • 

0,373 

’  -0.810 

-0*090 

-0.020 
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0.210 

B9o.5 

3447 
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9  0.023 

“0,00 
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0,06 
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0.208 
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3467 
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0.335. 

0^219 
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0 

0. 
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0.212 

8*0.5 

0 

0, 
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0,09 

0,474 

-0.778 

-0,034 

0,027 

0.032 

0.207 
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0 

0. 
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-0 , 02 
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4,41 

-0.0 

5,42 
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0 

0, 
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4,46 

-0,0 
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0,13 
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-0.031 
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0.208 
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0 

0, 
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Table  1.  Continued 


RUN  TRAJ 

Moc  P, 

T 

q 

*,  P« 

R*~  Tdp  sh 

A 

h 

St 

DATE 

TIME  CON  SET 

ZERO  SET  TRANSONIC  4T 
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,3  541 

,  4  377 

B  9  749,3 

4  73  . 
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55  s  46 

151/  2 

146/ 

2  TEST 

TC 

-635 

STORE 

Wt 

A 

P, 

f2'P3 

X  „ 

AX  AX  Y 

1 

XX 

'xy 

. 

C„ 

c 

r 

mKt83 

eg 

in, eg 

n,cg 

eg 

eg 

XZ  'YY  YZ 

'zz 

rp 

m  q 

Lnr 

985..  1 

.069 

1.16/ 

1.167  4 

042 

0  ,  000  a 

,000  U. 

000 

o.uoo 

4.6 

u .  o 

3fU  106 

.6  0.8 

106,6 

-2.9 

-85.9 

- 

85.9 

A/C 

a  P 

Nz 

Y 

_^A/C  *P 

1 

Y  !R 

CONFIG 

.WING 

MOTION 

NOROLL 

POST 

COEF  THRUST  EJECT  '  XFE 

a'xae 

OJm 

?ei 

ZE2 

F-4E  1 

98  U.IO 

1.0 

0.0 

0.0  -1,00  0  , 

DO  0,0 

3 

FUSPL 

6 

i~ 

u 

O 

0  2 

3.46 

1.50 

"u'o  0 

38 

0.38 

FULL 

SCALE  VELOCITIES  AND 

accelerations 

PN  t 

VX 

VY 

vz 

Vr 

U 

v 

w 

P 

r 

5 

V 

W 

P 

9 

r 

5  0.003 

932,4 

u  .u 

15.9 

932,5 

0,  u 

u.u 

U  ,0 

0,00 

U,00 

o.oo 

-7.1 

-2.8 
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-0.30 
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TO  .0 
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0 .0 
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0  o  26 
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-12.2 

0.1 
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0.1 
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12  0.05; 
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-20  *  0 
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-27.40 

-0.14 

13  0.06’ 

930,7 

0.3 
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-0.8 
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0,4« 

0,72 

-0.01 
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-8.35 
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0.4 
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“0.01 
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932.2 

-1.6 

0.5 

16.  i 

o»ai 

0,30 

-0.01 

-13,6 

11.9 

16.4 

6,47 

-8,33 

0.04 

19  0,12., 

928,1 

1.5 
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932.1 

-!•  / 

0 ,6 

16.3 

0,88 

0.22 

-0.01  ■ 

“12,4 

13,2 
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7.09 

-8.24 

0.10 
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927  t  9 
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88.1 

932,1 

*1.9 

0  o  7 

16,4 

0,95. 

0,14 

-0.01 
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16,3 

6,74 

-8.19 

0.17 
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2.1 
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Table  1.  Concluded 


U) 


RUN  TRAJ  p,  Tf  q_  p„  T_  ReM  Tpp  SH  A  li  Si  DATE  TIME  CON  SET  ZERO  SET  TRANSONIC  4T 

151  1312  0,B«9  1200,3  541.4  317, V  749,3  473,2  2.4  457.3  0,0012  0-050  5, OK  0,0010  11727/79  19l55:46  151/  2  146,-  2  TEST  TC-635 

ST0RE  W>  A  ri  f2'f3  X<3  AXn,.c  AX..c3  V  lxx  Ixt  ixz  lry  Iyz  lzz  Qp  C.,  C„, 

MKt63  »85i:  1,069  1,167  1.167  4,042  0,000  0.000  U,000  0,000  4,8  0,0  0.0  106,9  0,8  IBS.,6  -2.9  -B5.9  -85.9 

A/C  a  p  Nz  Y  <3A/C  lp  lY  lR  CONFIG  WING  MOTION  NOROLL  POST'  COEF  THRUST  EJECT  XFE  AXae  ZE)  ZE2 

F-4E  1.9«  0.00  1.0  0.0  0.0  -1,00  0.00  0,0  3  FUS?L  0  1  0  0  0  2  3,46  1,50  0.0  0,38  0.38 


PYLON  AXIS  POSITIONS  AND  orientations' 
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ie 
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19 
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"0.1 
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20 
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0,0 

1.6 

-0.1 

3.4 

-0,0 

21 
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-0 , 0 
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Table  2.  Standard  Aerodynamic  Grid  Tabulated  Summary  Data  Format 
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SURVEY  Mw 

p, 

T, 

<? 

P» 

Toc 

SH 

A 

DATE 

TIME 

CON 

SET  ZERO  SET  TRANSONIC 
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1203, i  535,8 
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.2 
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7  650,2  2, 
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a 

|8 

*p 

1, 
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Y 

eg 

Z 

eg 

F/S 
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0.0 
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1 
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MK 
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c„ 
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12 
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-0 
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5 
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0 

.029 

0,128 

237.1 

6 

166 
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7 
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0.01 

7,85 

-0 

,02 

1.1 

0 

*12 

*7,85 

0,063 

-0,079 
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-0.C2 

?o,oc 
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0 
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237.5 

14 
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24 
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-0,01 
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6 
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6 
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0 
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25 
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0.01 
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7 
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8 

,27 

*7,40 

0.662 

-U.862 

0,587 

-0.929 

0 

.046 
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237,5 

16 

166 

2 

-0,02 

-0.02 

"0.01 

-o.oo 

0 

,05 

-0.0 

8 

,05 

0,0  4 

0,565 

-U.585 

-0,010 

O.01Z 

0 

,02? 

0/124 

237.2 

17 

167 

3 

-0,02 

0.00 

-0.02 

o.ol 

3 

,98 

0.0 

3 

,98 

-0.01 

0,286 

-0,280 

-0 . 013 

0,064 

0 

,030 

0/120 

237.3 

18 

167 

4 

-0.03 

.0,01 

-0 . 02 

0,01 

1 

,99 

0.0 

1 

,99 

-0,01 

0,163 
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-o.oi3 

0,05/ 

0 

,032 

0/133 

238,2 

19 

-167 

5 
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0,03 
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.00 

0.0 

1 

,00 
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0.04/ 
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.030 
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20 
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0.01 
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0.01 
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.00 
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0/134 

230,2. 

_ 21- 

_ 147.  - 

7 

-0.00 

o..u: 

"0.02 

0,01 

-1 

■  00 

0.0 

-1 

,00 

•0,01 

-6.023 

U  ,  039 
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22 
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a 
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-2 

,00 
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-2 
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0 
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23 
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9 
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0.01 

-4 
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0.0 

-4 

,00 
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-0.214 
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-0,024 

0 

.034 
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24 

167 

10 

-O.ol 
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0.00 

-7 
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-7 
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25 
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29 
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15 

-0,02 

0.00 
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'0,00 
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-0.0 

-28 

,06 

-0,00 

-2,499 
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-0,151 
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0 
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237,2 

30 
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-o  ,  01 

:o.oi 

-0,0  3. 
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,05 

0.0 

-32 

,05 

0,02 

-3,138 

2,532 
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0,727 

0 

,042 
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31 
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17 

-0.C3 
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*36 

,11 

0.0 

-36 

oil 

0,05 

-4,09/ 

3.115 
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Table  3.  Standard  Flow-Field  Grid  Tabulated  Summary  Data  Format 


RUN  SURVEY  M_  P,  T,  ^  T_  V„  R«m  Top  SH  k  DATE  TIME  CON  SET  ZERO  SET  TRANSONIC  4T 

979  131  0.952  1740.5  96.6  616.4  971.6  470.9  1012.7  3.5  10.5  0.0017  0.050  11/21/79  16/31/27  877/  3  019/  2  TEST  TC-623 

A/C  a  y3  Ip  lY  lR  CONFIG  WING 

F-4C  17.00  0.00  -1.00  0.00  0.0  35  LEFT 


SUMMARY  1 


reference  axis  probe  pressures  flow  angles 


RUN/PN 

XREF 

Y 

1  REF 

ZREF 

p..l 

p.,2 

pO 

P.,4 

PP,5 

p5 

p5,P 

Ml 

Cp, 

f 

CP,F 

a 

NDX 

879003 

-14.75 

-0.02 

1.17 

-0.02 

-0.02  -0.1 

1280.1 

1309.5 

1178,7 

1169.0 

1724.0 

0..716 

0.7231 

0.82 

0.194 

5.4 

-0.268 

-6.5 

51 

879004 

-14.24 

-0.03 

1 .18 

-0.02 

-0.02  -0.1 

1283.2 

1311.1 

1182.2 

1172.0 

1723.8 

0.718 

0.7247 

0.02 

0.194 

5.4 

-0.267 

-6.5 

52 

87900? 

-13.75 

-0.03 

1.17 

-0.03 

-0.02  0.0 

1288.8 

1318.6 

1188.1 

1177.4 

1724.1 

0.721 

0.7283 

0.B1 

0.196 

5.5 

-0.275 

-6.7 

53 

879006 

-13.25 

-0.03 

1.17 

-0.02 

-0.02  -0.1 

1290.9 

1322.1 

1190.4 

1178.9 

1724.0 

0.722 

0.7298 

0.81 

0.196 

5.5 

-0.280 

-6.8 

54 

879007 

-12,74 

-0.03 

1.18 

-0.03 

-0.02  -0.1 

1292.7 

1327.0 

1192.4 

1180.5 

1724.3 

0.724 

0.7314 

0.80 

0.197 

5.5 

-0.288 

-7.0 

55 

879008 

-12.24 

-0.03 

1.18 

-0.03 

-0.02  0.0 

1297.0 

1336.3 

1196.8 

1184.4 

1724.0 

0.727 

0.7352 

0.79 

0.200 

5.5 

-0.304 

-7.4 

56 

879009 

-11.75 

-0.03 

1.17 

-0.02 

-0.02  -0.1 

1302.7 

1342.6 

1202.1 

1189.0 

1724.0 

0,730 

0.7304 

0.79 

0.203 

5.6 

-0.310 

-7.6 

57 

879011 

-11.24 

-0.03 

1.18 

-0.03 

-0.02  -0.0 

1305.7 

1349.3 

1204,5 

1191.1 

1724.4 

0.732 

0.7407 

0.78 

0.206 

5.7 

-0,322 

-7.9 

58 

879012 

-10.74 

-0.03 

1.18 

-0.03 

-0,03  0.0 

1308.6 

1355.0 

1205.7 

1192.0 

1724.3 

0.734 

0.7427 

0.77 

0.21 1 

5.8 

-0.335 

-0.2 

59 

079013 

-10.25 

-0.02 

1.17 

-0.01 

-0.02  -0.0 

1311.4 

1358.0 

1205.9 

1192.3 

1724.3 

0.735 

0.7440 

0.77 

0.218 

6.0 

-0,342 

-8.4 

60 

079014 

-9.74 

-0.03 

1.18 

-0.02 

-0.03  0.0 

1315.9 

1362.5 

1206,4 

1193.5 

1724.4 

0.736 

0,7460 

0.77 

0.228 

6.3 

-0.351 

-8.7 

61 

079015 

-9.25 

-0.03 

1.17 

-0.02 

-0.02  -0.1 

1320.5 

1366.4 

1206.4 

1194.0 

1724.4 

0.738 

0.7478 

0.76 

0.239 

6.5 

-0.361 

-8.9 

62 

879016 

-8.74 

-0.03 

1.19 

-0.02 

-0.03  0.0 

1324.6 

1368.1 

1203,0 

1191.7 

1724,6 

0.737 

0.7486 

0.76 

0,256 

7.0 

-0.371 

-9.2 

63 

879017 

-8.25 

-0.03 

1.17 

-0.02 

-0,02  -0.1 

1329.4 

1369. 8 

1201 ,0 

1190.5 

1724.3 

0.738 

0.7499 

0,76 

0.271 

7.4 

-0.379 

-9.4 

64 

879010 

-7.75 

-0.02 

1.17 

-0.02 

-0.02  -0.1 

1334.1 

1370.7 

1197,7 

1187.9 

1724.4 

0.738 

0,7506 

0.76 

0,289 

7.8 

-0.387 

-9.6 

65 

879019 

-7.25 

-0.02 

1.17 

-0.01 

-0.02  -0,0 

1338.3 

1370.8 

1193.1 

1184.4 

1724.5 

0,737 

0,7510 

0.75 

0.308 

0.3 

-0.396 

-9.9 

66 

879020 

-6.75 

-0,03 

1.17 

-0.02 

-0,02  -0.1 

1338,9 

1367.5 

1184,6 

1176,7 

1724.7 

0,735 

0.7492 

0.76 

0.325 

8.7 

-0.402 

-10.0 

67 

879022 

-6.24 

-0,03 

1.18 

-0.02 

-0.02  -0.1 

1341.1 

1363.0 

1178.8 

1171.8 

1724.7 

0.733 

0.7480 

0.76 

0.340 

9.1 

-0.401 

-10.0 

68 

879023 

-5.75 

-0.03 

1.17 

-0.03 

-0.02  -0.1 

1340.8 

1359.3 

1172,2 

1165.4 

1724.4 

0.730 

0.7464 

0.77 

0.351 

9.4 

-0.404 

-10.1 

69 

879024 

-5.25 

-0.02 

1.17 

-0.02 

-0.02  -0.1 

1343.2 

1354.8 

1169.2 

1162.3 

1724.6 

0.729 

0.7455 

0.77 

0.361 

9.7 

-0.399 

-10.0 

70 

079025 

-4.74 

-0.02 

1.18 

-0.02 

-0.03  0.0 

1347.5 

1.353.1 

1169.8 

1162.6 

1724.8 

0,730 

0.7461 

0.77 

0,370 

9.9 

-0,396 

-9.9 

71 

879026 

-4.24 

-0.02 

1.18 

-0.01 

-0.03  0.0 

1353.3 

1349.7 

1 170 i 9 

1163.4 

1724.8 

0,730 

0.7468 

0.76 

0,380 

10.2 

-0.389 

-9,7 

72 

879027 

-3.74 

-0.03 

1.18 

-0.02 

-0,02  -0.1 

1356,0 

1344.0 

1167,2 

1159.7 

1724.7 

0,729 

0.7458 

0,77 

0.392 

10.5 

-0.383 

-9.6 

73 

879028 

-3.25 

-0.03 

1.18 

-0.02 

-0.03  0.0 

1363.8 

1339.9 

1170.1 

1162.3 

1724.7 

0.730 

0.7471 

0.76 

0.404 

10.8 

-0.371 

-9.3 

74 

879029 

-2.75 

-0.02 

1.18 

-0.01 

-0.03  0.0 

1372.4 

1335.5 

1173,3 

1164,6 

1724,7 

0.731 

0.7485 

0.76 

0,418 

11.2 

-0.359 

-9.0 

75 

879030 

-2.25 

-0.02 

1.17 

-0.02 

-0.02  -0.1 

1382.2 

1328.6 

1172.5 

1163,1 

1724.5 

0.732 

0.7494 

0.76 

0,442 

11.8 

-0.349 

-8.8 

76 

879031 

-1.75 

-0.03 

1.18 

-0.02 

-0.02  -0.1 

1387.2 

1316.6 

1164,5 

1155,5 

1724,6 

0,728 

0.7472 

0.76 

0,465 

12.4 

-0.337 

-8.5 

77 

879032 

-1.25 

-0.03 

1.18 

-0.03 

-0.02  -0.1 

1390.3 

1300.4 

1152,2 

1143.8 

1724.4 

0.723 

0.7434 

0.77 

0,490 

13.0 

-0.322 

-8.1 

70 

879033 

-0.74 

-0.03 

1.18 

-0,02 

-0.03  -0.0 

1394.4 

1287.5 

1142.7 

1135,1 

1724.5 

0.719 

0.7408 

0,78 

0.513 

13.6 

-0.310 

-7,8 

79 

879034 

-0.25 

-0.02 

1.18 

-0.02 

-0.02  -0.0 

1394.9 

1273.3 

1130.4 

1124,7 

1724.6 

0.714 

0.7369 

0.79 

0.531 

14.1 

-0.298 

-7.5 

80 
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Table  3.  Concluded 


RUN 

SURVEY 

Mao  pt 

T, 

p«>  T, 

Re^ 

TOP‘ 

SH 

A 

DATE 

TIME 

CON  SET 

ZERO  SET 

TRANSONIC  4T 

925 

136  0 

.603  2219. 

1  96.5  441.5 

1736.1  510.5  672.7 

3.5 

30.4 

0.0034 

0.050 

11/21/79 

19/21/49 

921/  12 

912/  1 

TEST  TC-623 

A/C 

a 

ft 

'p  'y  'r 

CONFIG  WING 

F-4C 

5.00 

0.00  -1 

.00  0,00  0 , 

0  55  LEFT 

SUMMARY  2 

PYLON 

PN  Xp 

AXIS 

YP 

ZP 

“XY 

BODY  AXIS 

aXZ  aYZ 

FLOW 

vx 

ANGLES  . 

V 

v  X  Y 

AND  VELOCITIES 

VXZ  VY 

VYZ 

VZ 

PM> 

VL 

"t 

ml 

NDX 

RUN 

1 

-5.84 

0.33 

0 . 36 

-3.4 

1.0 

286.5 

650.1 

651.2 

650.2 

-38.7 

40.4' 

11.5 

2216.6 

418.0 

651.3 

3.5 

0.58 

101 

925 

2 

-5.55 

0.33 

0.36 

-3.5 

1.3 

290.4 

647.9 

649.1 

640,0 

-40.1 

42.8 

14.9 

2216.7 

415.8 

649.3 

3.8 

0.58 

102 

925 

3 

-5.25 

0.33 

0.35 

-3.7 

1  ;7 

295.1 

651.7 

653.1 

652.0 

-41.7 

46.0 

19.5 

2220.0 

420.8 

653.3 

4.0 

0.58 

103 

925 

4 

-4.95 

0.33 

0.36 

-3.8 

2.1 

298.6 

655.8 

657.3 

656.3 

-43.5 

49.5 

23. -7 

2215.9 

424.7 

657.7 

4.3 

0.59 

104 

925 

5 

-4.65 

0.33 

0.36 

-3.9 

2.3 

301.1 

660.0 

662.3 

661.3 

-44.7 

52.2 

27.0 

2215.7 

430.2 

662.8 

4.5 

0.59 

105 

925 

6 

-4.35 

0.33 

0.35 

-3.9 

2.5 

302.9 

663.3 

664.9 

664.0 

-45.3 

54.0 

29.4 

2214.9 

432.8 

665.5 

4.6 

0.60 

106 

925 

7 

-4.05 

0.33 

0.35 

-3.9 

2.5 

303.1 

671.2 

672.8 

671.9 

-.45.6 

54.4 

29.7 

2214.8 

441.3 

673.4 

4.6 

0.60 

107 

925 

8 

-3.75 

0.33 

0,36 

-3.8 

2.3 

301.4 

675.8 

677.3 

676.4 

-44.6 

52.2 

27.2 

2213.9 

445.8 

677.8 

4.4 

0.61 

108 

925 

9 

-3.45 

0.33 

0,35 

-3.6 

1.9 

297.2 

677.3 

678.7 

677.7 

-42.7 

48.0 

21.9 

2215.0 

447.4 

679.0 

4.0 

0.61 

109 

925 

10 

-3.15 

0.33 

0.35 

-3.4 

1.2 

289,9 

679.6 

680.0 

679.7 

-40.8 

43.4 

14.8 

2216.0 

449.5 

680.9 

3.6 

0.61 

110 

925 

11 

-2.85 

0.33 

0.36 

-3.2 

0,6 

280.3 

6B0. 4 

681,5 

680.5 

-38.3 

38.9 

6.9 

2217.3 

450.4 

681.5 

3.3 

0.61 

111 

925 

12 

-2.55 

0.33 

0.36 

-3.0 

0.0 

270.6 

677.9 

678.8 

677.9 

-35.5 

35.5 

0.4 

2214.6 

447.0 

678.8 

3.0 

0.61 

112 

925 

13 

-2.25 

0.34 

0.35 

-2.8 

-0.4 

261.1 

676.4 

677.2 

676.4 

-33.0 

33.4 

-5.2 

2215.3 

445.5 

677.2 

2.8 

0.61 

113 

925 

14 

-1.95 

0.13 

0.35 

-2.6 

•1.1 

246.9 

674.0 

675.5 

674.9 

-30.2 

32.9 

-12.9 

2217.0 

444.1 

675.6 

2.8 

0.61 

114 

925 

15 

-1.65 

0.34 

0,35 

-2.3 

-1.7 

233.2 

658.6 

659.1 

658.9 

-26.3 

32.9 

-19.7 

2218.7 

427.0 

659.4 

2.9 

0.59 

115 

925 

16 

-1.35 

0.33 

0.35 

-2.2 

-1.9 

229.0 

654.2 

654,6 

654.5 

-24.6 

32.7 

-21.4 

2212.4 

421.2 

655.0 

2.8 

0,59 

116 

925 

17 

-1.05 

0.34 

0.35 

-2.2 

-1.7 

232.7 

647.4 

647.9 

647.7 

-24.7 

31.0 

-18,8 

2213.2 

414.1 

648.2 

2.7 

0.50 

117 

925 

18 

-0.75 

0.313 

0.36 

-2.2 

-1.2 

241.4 

642.0 

643.3 

642.9 

-24.7 

28.1 

-13.5 

2215.4 

409.4 

643.4 

2.5 

0.57 

118 

925 

19 

-0.45 

0.34 

0.35 

-2.2 

-0.6 

254,9 

633.5 

634.0 

633.5 

-24.5 

25.4 

-6.6 

2213.3 

399.0 

634.0 

2.3 

0.57 

119 

925 

20 

-0.15 

0.33 

0,35 

-2.3 

0.2 

274.3 

634.7 

635.2 

634.7 

-25.2 

25.2 

1.9 

2210.5 

399.7 

635,2 

2.3 

0.57 

120 

925 
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APPENDIX  A 

CONSTANTS  ASSEMBLY  PROCEDURE,  WITH  STAGING  AND 
INITIALIZATION  EQUATIONS 

A-l.  CONSTANTS  ASSEMBLY  PROCEDURE 

To  store  a  unique  set  of  the  360  constants  required  for  every  proposed  trajectory  in  a 
typical  wind  tunnel  test  program  is  not  practical  for  many  reasons.  Therefore,  a  “constants 
merge”  program  was  developed  to  assemble  the  necessary  360  constant  inputs  into  a 
temporary  array,  using  a  categorized  permanent  storage  file.  The  merge  routine  is  activated 
each  time  a  constant  point  is  requested  and  updates  the  temporary  file  which  is  to  be  used  in 
trajectory  calculations. 

The  structure  and  operation  of  the  merge  program  are  outlined  as  follows.  A  disk  file  is 
partitioned  into  eight  different  groups  —  A,  B,  C,  D,  E,  F,  G,  and  H.  The  “A”  file  is 
designated  the  main  deck  group  which  allows  from  1  to  5  arrays  of  360  constants  each.  The 
constants  which  are  essentially  invariant  for  the  test  are  stored  in  the  “A”  file.  Normally, 
only  one  main  deck  is  necessary  unless  more  than  one  balance  is  to  be  used.  The  “B,”  “C,” 
and  “D”  files  each  allow  from  1  to  9  .arrays  of  60  constants  while  the  “E,”  “F,”  “G,”  and 
“H”  files  each  allow  from  1  to  99  arrays  of  60  constants.  These'  files  are  used  to  store 
constant  values  which  change  during  testing  (e.g.,  altitude,  ejector  forces).  These  “variable” 
constants  are  divided  into  compatible  groups  as  determined  by  test  requirements,  and  each 
group  is  assigned  a  numerical  location  in  one  of  the  “B”  through  “H”  files.  Using  the 
merge  program,  one  can  now  recall  and  assemble  the  constants  groups  through  their 
alphabetical  and  numerical  designations  to  provide  the  unique  set  of  inputs  for  specific 
trajectory  applications. 

Prior  to  starting  a  trajectory,  the  numerical  values  of  each  group  required  are  dialed  into 
the  respective  “A”  through  “H”  constants  boxes  on  the  data  acquisition  panel  (DAP). 
When  a  constant  point  is  requested,  the  program  scans  these  constant  boxes  to  define  which 
arrays  are  to  be  assembled  and  then  loads  the  appropriate  groups  into  the  temporary  file  in 
alphabetical  sequence.  After  the  staging  calculations  are  performed,  the  calculated 
parameters  are  appended  to  the  file,  which  is  then  ready  for  access  by  the  trajectory 
generation  job. 
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A-2.  STAGING  EQUATIONS 

Calculations  done  during  the  staging  process  are  listed  below: 
Model  Reference  Area  and  Lengths 

Am  ,  A«2)(,.a)2 

f,  „,  .  f,  (12)  (A)  «kA) 

f2,m  '  l'2(12HA)(tA) 

V  -  (' 

1  3  ,m  “  1 2,m 

Model  Moment  Transfer  Distances  (See  Fig.  A-l) 

Xm  =  Xcg  (12)  (A)  (kA) 

•  Yni  =  Ycg  (12)  (A)  (kA> 

,Zni  =  Zcg(12)(A)(kA)  • 

Xm,t  =  XBF+Xm,ec-Xm  +  AXmc6(12)(A)(kA) 

Xn,t  =  XBF  +  Xn,ec  -Xm  +  AXn,cg  (12)  (A.)  (k^) 


CTS  Rig  Physical  Parameters  (See  Fig,  A-2) 

^l.R  =  |[(XP1TC  +  XBF  _Xm}//cos0sJ  ~3  j/cos'/'s 

^2,R  =  _(f  1 , R  sin  *As  +  YPITc) 
d p  =  (fj  j|  c os  +  3)  sin  0g  +  ^pjrp£ 

f3,R  =  (^,R  +  (Xm  cos0S"dR  sin0S>  cos^S  +  £2,R  sin^S 

Aircraft  Angle  of  Attack 

a  =  “Ct'pp 
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Wing  Flag 

If  (Ytp  >  RW),  Wing  =  Right  (1) 
If  (LW  <  YTp  <  RW),  Wing  =  Fusel  (2) 
If  (YTp<LW),  Wing  =  Left  (3) 

Store  Initial  Attitude  -  Launch  Trajectory 


l'R,I 

=  a  +  I  p  - 

-*s 

7r,I 

=  Iy-/3- 

-  i/js 

WR,I 

=  !r 

[**} 

=  TME  (u 

R,I*  ^R.I’  wR,I’  !'* 

[bg] 

=  TME  (i/i 

'S’  0S'  °>  2) 

[«■] 

-  M  j 

M 

’ll 

=  sin-1  [ 

CC  (1,2)] 

VY 

=  sin-1  [ 

— CC  (l,3)/cOS  7]l  ] 

cn, 

=  tan-1  [ 

-CC  (3,2)/CC  (2,2)] 

Store  Initial  Attitude  -  With  Induced  Angle  Corrections 

When  uo,  v Q,  wo  4  0  : 

h'  =  6.356766(103)  h /  [2.0855532  (107)  +  h] 
Ta  =  389.97 

If  (h  '<  11),  Ta  =  518.67  -  11.7  (h') 

♦See  Appendix  M  for  matrix  definitions. 
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UA  =  49.021  (Mg) 


v\  =  v\ 


m  =  »7i 


CO  j  =  CO  ( 

If  (POST  #  0),  vj  =  ^I)0,  J?|  =  T/i  0,  &>i  =  wI>0 
=  TME  (fj,  77 j,  6>j,  1) 


XI,o 

U  ^ 

r  1' 

O 

Ko 

=  f  BB 1 

V 

L  J 

0 

Z, 

w 

1,0 

0 

I't  =i/]+tan  1  [zIi0/(XIj0  +  UA)] 
r?[  =  i/i-tan-^Yj  o/[(XI  o  +  UA)  cos  ul  - \0  sin  u] 

CO  t '  =  6)j 


Store  Initial  Position  -  Unaccelerated  Flight  (Nz  =  1) 


X,  =  XT 
t  1,0 


Y,  =  Y, 
t  1,0 


=  Z 


1,0 


Store  Initial  Position  -  Accelerated  Flight  (Nz  ^  1) 

az  =  (1  _Nz)  32-174 


-a  7  /Ua 
Z  ,p  A 
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R 


P 


Z.p 


q 


2 

P 


0p  =  q p(t 0)  (57.2958) 

XC  =  XI,o+  RPsin01>+  lJAlo 
YC  =  YI,o 

ZC  =  Zl.o-V1  _cosV 


_AA]  =  TME  (0,-0  ,  0,  2) 


xt 

I-  -1  - 

XG 

Y, 

<r 

< 

II 

YC 

zt 

zc 

A-3.  INITIALIZATION  PROCEDURE 


Integrators  are  assigned  and  initialized  as  part  of  the  initialization  routine.  Integrator 
outputs  (designated  PI  through  P50)  are  assigned  initial  values  as  indicated  below. 

Initial  Inertiai-to-Body-Axis  Matrix 

^R,I  =  a  +  Ip  -fls 


'iR.I  =  ly  ~  P  ~^S 
WR,I  =  JR 

[AA  j  =  TME  (vR>1,  p  <uR>p  1) 

[BB]  =  TME  (<frs,  9S,  0,  2) 

[tobodyJ  =  [bb]  [aa] 

If  (POST  4  0),  [tOBODYo]  =  TME  (^  D,  Vl  o,  col  o,  1) 
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Initial  Values 


Parameters 


PI  =  U0 

P2  =  vo 
P3  =  vvc 

P4  =  X  j 

P5  =  Y 
P6  =  Z 


1,0 


1,0 


P7  «  Po 

P8  =  % 
P9  =  r„ 


P10  PI 3  P16 
Pll  P14  P17 
PI 2  PI 5  PIS 


P19  =  0 
P20  =  0 
P21  =  0 


Body  — axis  linear  velocities 


Inertial —axis  positions 


Body —axis  angular  velocities 


[TOBODyJ  Ine  rtial  —to  —body -axis  —matrix  direction  cosines 


Rotation  center  (hook)  velocities  (special  staged 
release  applications) 


Program  parameters,  flags,  and  counters  which  are  set  on  the  first-pass  loop  include  data 
cycle  time  (At),  number  of  integrations  per  data  cycle  (K;nt),  integration  time  (<5t  =  At/Kjm), 
initial  trajectory  time  (tQ),  integrator  pass  counter  (NPASS  =  0),  data  cycle  pass  counter 
(PASS  =  0),  extrapolator  pass  counter  (NXP  =  0),  number  of  extrapolators  (NX),  ejector 
control  flag  (EJECT),  staged  separation  control  flag  (MOTION)  and  data  cycle  time 
increase  control  flag  (STEP). 


Once-only  calculations  are  listed  below. 

Conditions  at  Simulated  Altitude* 

h  '  =  6. 356766(1 03)  h/[  2.0855532  (107)  +  h] 

If  (h '  <  11).  Ta  =  518.67-  11.7  (h') 

pA  =  2116.216  [l  -0.0225577 (h')]5,2559 


The  pressure  and  temperature  at  altitude  equations  are  valid  from  sea  level  to  65,800  ft.  Reference:  1975 
U.S.  Standard  Atmosphere  from  NASA-TR-R-459  (May  1976). 
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If  (h  '  >  11),  Ta  =  389.97 

pA  =  472.677  [e-°-157688  (h'-11)] 


UA  =  49.021  Mb>/T7 


P\  =  1-4PaMBUA 

Weight  Components,  Mass 


Wy  =  Wt  sin  y 


R  y  =  Wt  cosy  sin  ^  A/C 
wz  =  Wt  cosy  cos  cf,  A/c 
m  =  Wt/32.174 


Initial  Store  Attitude  (Pitch,  Yaw,  Roll  Sequence) 


sin-  1 

[tOBODYo  (1,2)] 

sin- ^  | 

[-TQB0DYo(1,3)/cos7o] 

tan-  *  j 

-T0B0DYo  (3,2)/T0B0DYo  i 

2,2)] 

Initial  Store  Angular  Rates  for  Accelerated 
Flight:  Launch  Trajectory  Only 


Po,p 

0 

^o,p 

=  [tobodyo] 

qp 

r 

o 

o,p 

P7  =  P7+Po,p 

P8  =  p8  +  q0,p 


P9  =  P9  +  r 


o,p 
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Inertia  Matrix 


Inverse  Inertia  Matrix 


[1] 


[I] 


•xx 

-IXY 

-!xz 

X 

-< 

iyy 

_IYZ 

-*xz 

_IYZ 

}ZZ 

- 1  = 

inverse  of 

[I] 
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CJN 

on 


MOMENT 

REFERENCE 

CENTER 


CG 


STORE 

CG 


REAR  VIEW 


SIDE  VIEW  NOTE  :  AXm^AXn.cg  ARE  USED  TO  TRANSFER 

THE  Cm  AND  Cn  MOMENT  REFERENCES 
FROM  THE  CG  POSITION.  HOWEVER, 
MOTION  IS  STILL  CALCULATED  ABOUT 
THE  TRUE  CG  POSITION  AND  THE  OFFSET 
COEFFICIENT  DATA  ARE  USED  AS  THOUGH 
OBTAINED  ABOUT  THE  TRUE  CG  POSITION. 


Figure  A-1.  Store/balance  physical  definition. 
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CG  -  0S 


BAL  (L 


CTS 

CTS  PITCH 
YAW  CENTER 
CENTER 


SIDE  VIEW 


BAL  I. 


M 


CTS 

CTS  PITCH 
YAW  CENTER 
CENTER  [ 


TOP  VIEW 


DIMENSIONS  IN  INCHES 

NOTE:  BALANCE  <£.  AT  ZERO 
PITCH  AND  YAW 


Figure  A-2.  Captive  trajectory  support  (CTS)  mechanism 
physical  definition. 
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Angular  Sting  Deflections 

A v  =  [Ka,NF>T,g  +  Ka,mMm,g]C0SwR  ~[K^,YFY,g+  K^>nMn,g]  sin&JR 
A7  =  [K^,YFY,g  +  K(A,nMn,g  ]C0SwR  +  [Ka,NFN,g+  Ka,mMm,g]  sinajR 

Aw  =  K^,fMelg 

Linear  Sting  Deflections 


axb  =  0 

AY  B  =  KY,Y 

[FY,g  cos  ^R  "+  FN , g  s‘n  WR 

l+KY,n| 

[Mn,g  C0SWR  +  Mm,g  sinA>R] 

azb  =  KZ,N  [ 

_FN,g  cos  WR  _  F  Y,g  sinwR 

|  +  KZ,m  | 

[Mm,g'cos<uR-Mn,gsin<uR] 

[AA]  =  TME  („t,  Vl,  0,  1) 


axl 

axb 

ayl 

II 

1 - 1 

> 

> 

1 - 1 

ayb 

azl 

AZb 

Note:  Typical  examples  showing  how  angular  and  linear  sting  deflection  constants  are 
obtained  from  static  loadings  are  presented  in  Fig.  B-l. 

Static  Tare  Corrections 

[AA]  =  TME  (v(,  t]v  cov  1) 


AAA 

0 

BBB 

=  [AA] 

0 

CCC 

1 
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If  (SRj 

=  0  and  <5RS  =  0)  CCC  =  CCC  -  1 

F’  N.st 

=  -CCC(WN) 

-  # 

fy  , 

I  ,st 

=  BBB  (\VY) 

' 

^  A ,  s  t 

=  -AAA  (WA) 

- 

Mp.Bl 

=  CCC  (Y)  WN  -BBB  (Z)  Wy 

Mm.Ht 

=  AAA  (Z)WA-  CCC  (XJWN 

M  , 

n,st 

=  BBB  (Xn)  Wy -AAA  (Y)  WA 

Net  Loads 

fn  = 

FN,g-FN,st' 

FY  = 

FV,g  ~  FY,st  ' 

fa.,  = 

■ F  A ,  g  ~  F  A ,  s  t 

mf  = 

MP,g“MP,st  *  Ym(FV  +  Zm(FY) 

M  = 

rn 

Mm.g'-M»..t-Xn..|(FN)+Zm(FA.t) 

M  = 

n 

Mn,rM^-XM(FY)-Ym(FA,t) 

- 

Aerodynamic  Coefficients 

CN  =  FN/CIc<,Am 

CY  =  FY/cIooAm 

C  A,t  =  FA,Aoo  A  m 

CP  =  Mp/q«  Am  h,m 

Cm  =  Am  *l,n. 

“ 

Cn  =  M„/qM  Am 

- 
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Identify  Extrapolation  Parameters 


i  ,X 


=  C, 


=  C. 


'A,  t, 


=  c: 


A  ,t 


Cp 


Cf 


c  =  c 

n,x  n 
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CALIB  BODY  SLOTS 
O  SLOT  I 
□  SLOT  5 


.  FN»  # 

POSITIVE  IN  SAME  DIRECTION  AS  +v 


BALANCE 


SLOT  ELECTRICAL  SLOT 

I  CENTER  5 


DISTANCE  FROM  FACE 
OF  CALIB  BODY,  in. 


a.  Angular  (pitch  plane) 

Figure  B-1,  Typical  sting  deflections  constants  definition. 


W 


AYa/ib 
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CALIB  BODY  LOAD  SLOTS 
O  SLOT  I 
O  S  LOT  5 


AYb  POSITIVE  IN  SAME  DIRECTION  AS  +Fy 


SLOT 

I 


BALANCE  SLOT- 

ELECTRICAL  5 

CENTER  I 


b.  Linear  (side  plane) 
Figure  B-1.  Concluded. 
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APPENDIX  C 

- 

CONVERSION  MODULE  EQUATIONS 

Pickup  Integration  Results 

u  =  PI 

\ 

1 

V  =  P2 

>  Body —axis  linear  velocities 

w  =  P3 

\ 

Xj  =  P4 

Y,  =  P5 

Zj  =  P6 

Inertial —axis  positions 

p  =  P7 

I 

q  =  P8 

Body —axis  angular  velocities 

r  =  P9 

V 

[tobody]  - 

P10  P13  P16 

Pll  P14  P17 

PI  2  P15  P18 

Inertial  —to  —body  —  axis 

direction  cosine  matrix 

* 

u  R  =  Pl9 

=  P20 

w  B  =  P21 

^  Rotation  center  (hook)  velocities 

^  (staged  release  applications) 

- 

Euler  Orientation  (see  Fig.  C-l) 

0]  =  sin  1 

[  -TOBODY  (1,3)] 

i//j  =  sin-  1 

[  TOBODY  (l,2)/cos  0, J 

0  [  =  tan-1 

[  TOBODY  (2,3)  'TOBODY  (3,3)] 

Pitch,  Yaw,  Roll  Orientation  (see  Fig.  C-2) 

- 

r]  j  =  sin  1 

[  TOBODY  (1,2)] 

vj  =  sin  ' 

[  -TOBODY  (1,3)/ cos  77  j  ] 

* 

Cl>  |  =  tan  * 

[  -TOBODY  (3,2)  'TOBODY  (2,2)] 

- 
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If  (NOROLL  *  0),  call  nonrolling  capacity  routine  (see  Appendix  L.) 

Inertial-Axis  Velocities 


*1 

u 

Y, 

=  [tobody] 

V 

h 

w 

Weight  Components,  Total  Store  Velocity,  and 
Dynamic  Pressure  at  Altitude 


__ 

wx 

W  ' 

v'x 

wY 

=  [tobody] 

WY 

Wz 

WZ 

Y 

Ur  =  [(Ua  +  XjIMy^MZj)2]  2 
q  A  =  °-5PAlJR2 

Total  Velocity  Components  (Body  Axis) 


UA,X 

UA 

UA,Y 

=  [tobody] 

0 

lJA,z 

0 

V  y  =  U  a  y  +  v 

VZ  =  UA,Z  +  w 

Store  Angle  of  Attack  and  Sideslip  Angle 

ag  =  tan-  1  (V^/V^) 
|Ss  =  sin-  1  (VY/UR) 
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Lanyard  Length 


XL,o 

-X1 

Y  L.o 

=  [d'OBODY  ] 

“Y1 

Z  L.o 

-Z1 

dp  =  qp  (t)  57.2958 


X  L.S 

-X 

o 

YL,S 

=  [  tobody] , 

-Y 

o 

ZL,S 

-Z 

o 

Define  reference  point  on  aircraft  (initial  pass  only). 


Define  reference  point  on  store. 


AX 


L.S 


XL,S  +  XI  ~XL,o 


AY  L.S  ~  Y  L.S  +  X  I  ~  YL,o 


JL,S 


J  L.S  +  ^\~  ^L.o 


For  dp  =  0. 


AX 


AY 


AZ 


L.S 


L.S 


L.S 


XL,S  +  XI  ~  XL,o  +Rpsin^p+- 

yl,s-yi-yl,o 
ZL,S+ZI_ZL,o  -Rp(1-cos^p) 


•uAt  | 


For  0p  ?  0. 


1; 

ZL  c  =  (Ax l  g)2  +  (AYl  s>2  +  (^ZL  S)2  J  2  Lanyard  length. 
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-J 

on 


a).  FIRST  ROTATION, YAW  ANGLE 
(ROTATE  ABOUT  2j  AXIS) 


b)  SECOND  ROTATION,  PITCH  ANGLE(0x )  c)  THIRD  ROTATION.  ROLL  ANGLE  [<frx  ) 
(ROTATE  ABOUT  Yj  AXIS)  (ROTATE  ABOUT  X0  AXIS) 


Figure  C-1.  Graphic  illustration  of  a  yaw,  pitch,  roll  (Euler) 
orientation  sequence. 


AEDC-TR-79 


Xj  -  Zj  PLANE 


Zb  2i 

a)  FIRST  ROTATION, PITCH  ANGLE! Vj) 
(ROTATE  ABOUT  Yx  AXIS) 


Figure  C-2.  Graphic  illustration  of  a  pitch,  yaw, 
roll  orientation  sequence. 
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APPENDIX  D 

OFFSET  COEFFICIENT  MODULE  EQUATIONS 

The  flow  diagram  for  the  standard  offset  coefficient  module  is  shown  in  Fig.  D-l.  The 
equations  are  listed  below. 

Initial  Values,  Control  Flag  (Section  A) 

tDEL  =  0  cAo  =  (cAt0)o  cEi0=  (c£;0)o 

COEF  =  COEFI  CY  0  =  (CY  o)u  Cm>0  =  (Cm 
CN.o  =  (CNf0)o  Cno  = 

Ramp  Axial-Force  Equation  (Section  B) 

dC, 

CA,o  =  “dt  '  (t_tDEL)+  (C  A,o)q 
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Figure  D-1.  Flow  diagram  of  the  standard  offset  coefficient  module. 
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APPENDIX  E 

TOTAL  COEFFICIENT  MODULE  EQUATIONS 
Total  Coefficient  Calculations 

CA,t,T  =  CA,t,x  +  CA,o 

Cy,T  =  ^Y,x  +  ^Y,o  ' 

CN,T  =  CN,x  +  CN,o 

Cf,T  =  Cf\x  +  Cf,o  +  [pV2Ur] 

p 

C  T  =  C  +C  +C  r  qf\/2U„  1 

m ,  I  rn ,  x  m,.o  mq  L  *  V  RJ 

CT=C  +C  +C  f  rf„/2Un  1 

n,  I  n,x  n,o  nr  L  *  H  J 

If  (NOROLL,  =  3)  Cf  T  =  0 
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APPENDIX  F 

THRUST  MODULE  EQUATIONS 


The  flow  diagram  for  the  standard  thrust  module  is  shown  in  Fig.  F-l,  and  a  graphic 
description  is  presented  in  Fig.  F-2.  The  equations  for  calculating  thrust  force  and  moment 
contributions  are  as  follows: 


Initialization,  Control  Flags  (Section  A) 


THRUST  = 

T1IRSTI 

FT,X  = 

=  0 

^T,X 

tDEL  = 

0 

fT)Y  = 

=  0 

lVl^p  y 

tt.Cl  = 

0 

FT,Z  = 

=  0 

^T,  Z 

lt,C2  =  0 
Zt,C  =  0 

First  Polynomial  for  FTiX  (Section  B) 

FT.X  =  |  an  (tt)n 

Second  Polynomial  for  FTiX  (Section  C) 

'  ft,x  =  tw 


Thrust  Moment  Contributions  (from  Jet  Damping)  (Section  D) 

MT,X  =  Cjdf  (P)  ft,x 
mt  y  =  Cjd  (q)  ft,x 
MT  z  =  Cjd  W  F T,X 

J  n 
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THRUST  OPTIONS 

0  NO  THRUST 

1  NO  DELAY 

2  TIME  DELAY 

3  LANYARD  DELAY,  THEN  TIME  DELAY 


Figure  F-1.  Flow  diagram  of  the  thrust  module. 


EXIT 


8 


CONTROL  PARAMETER 


THRUST 

0  NO  THRUST  FORCES 

[  NO  DELAY  (TIME  OR  LANYARD  LENGTH) 

2  TIME  DELAY  ONLY 

3  LANYARD  DELAY,  THEN  TIME  DELAY 


Figure  F-2.  Graphic  description  of  thrust  force. 


APPENDIX  G 

EJECTOR  MODULE  EQUATIONS 


A  E  DC-TR-79-1 


The  flow  diagram  for  the  standard  ejector  module  is  shown  in  Fig.  G-l,  and  a  graphic 
description  is  presented  in  Fig.  G-2.  The  equations  for  calculating  forward  and  aft  ejector 
force  and  moment  contributions  are  as  follows: 


Initialization  Equations,  Control  Flags  (Section  A) 

El  FLAG  =0  X 


E2  FLAG  =  0 


X 


El 


E2 


'1C 


2C 


0 

0 


X  n  g  X  F  E 
XEl  ~  AXAE 


Set  Ejector  Forces,  Moments  to  Zero  (Section  B) 


fei  =  0 
^  E2  =  0 


E.X 


=  0 


fe,y  =  0 


me,x  =  0 


M 


E.Y 


^  E.Z  =  0 


ME,Z  =  0 


Stroke  Length  Equations  (Section  C) 

6p  =  qp  (t)  57.2958  - 
[AA]  =  TME  (fo  +  6p,  r]0,  0,  1) 


[BB]  = 

[tOBODy]  [AA]' 

Ad  = 

sin  *  | 

[-BB  (1,  3)] 

Ai p  = 

sin-  1  | 

[BB  (1,  2)/cos 

a@] 

Z1E  = 

[ZI  -RP  {1  -cos  V 

-  Xg  j  sin  Ad  J 

COS  (1] 

-  j^Y  |  +  X  E  j  sin  Aip  J 

sindp  +  ojm) 

Z2E  = 

[  Z1  “Rp  (]  -cos6»p) 

—  X  £2  sin  A 9 

jcos  (I 

-  [x  1  +  XE2  sin  A^]  sin  (Ip  +  com) 
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Forward  Ejector,  First  Polynomial 

5 


l'  El  =  £  Cn  (Z,E)n 

(Eject  <  2) 

fei  «  !c„Wn 

(Eject  =  3) 

Forward  Ejector,  Second  Polynomial 

FE1  =  |dn(ZlE)n 

(Eject  <  2) 

fei  = 

(Eject  =  3) 

Aft  Ejector,  First  Polynomial 

FE2  =  ?  en  ^Z 2E^n 

(Eject  <  2) 

FE2  =  |  en  Wn 

(Eject  =  3) 

Aft  Ejector,  Second  Polynomial 

FE2  ='iMZ2E)n 

(Eject  <  2) 

FE2'=!f„wn 

(Eject-=  3) 

Component  Resolution  Equations  (Section  D) 

E,X 

= 

0 

E,Y 

= 

~(FEl 

+  FE2)  sin  a>m 

E,Z 

= 

(fei  + 

F £2)  coswm 

^E,X 

= 

0 

lE,Y 

= 

-[fe, 

XE  1  +  FE2  XE2]  coswm 

E,Z 

= 

_[fei 

XE 1  +  FE2  XE2]sin 
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EJECT 

0  NO  EJECTOR  ACTION 

1  EJECTOR  FORCES  ■  f{t)  .CUTOFF  .  f(t) 

2  EJECTOR  FORCES  *  f  (STROKE).  CUTOFF  >  f  (STROKE) 

3  EJECTOR  FORCES  ■  f  {  t  )  ,  CUTOFF  ■  f  (STROKE) 


Figure  G-l.  Flow  diagram  of  the  ejector  module. 


85 


AEDC-TR-79-1 


EJECT  CONTROL  PARAMETER 

0  NO  EJECTOR  FORCES 

1  EJECTOR  FORCES  a  CUTOFF  =  f (TIME) 

2  .  EJECTOR  FORCES  8  CUTOFF  s f (STROKE) 

3  EJECTOR  FORCES  =  f (TIME), CUTOFF  =  f(STROKE) 


FWD  AFT 

STORE  EJECTOR  EJECTOR 

NOSE  PISTON  PISTON 


XFE 

■‘-^AE-* 

V _ 

+  CN 


ZB 


INITIAL  PORTION  FINAL  PORTION 

CURVE  FITS  CURVE  FITS 


Figure  G-2=  Graphic  description  of  ejector  forces. 
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APPENDIX  H 

FULL-SCALE  FORCE  AND  MOMENT  MODULE  EQUATIONS 

For  a  free-falling  store,  components  of  the  full-scale  forces  and  moments  acting  through 
or  about  the  store  eg  are  described  as  follows: 

Forces 

FX  =  A  CA,t,T  +  FT,X  +  E.X 

f  Y  =  +  ‘l A  A  FY,T  +  FT,Y  +  FE,Y 

fZ  =  ^Z  ~  ^A  A  CN,T  +  FT,Z  +  FE,Z 

Moments 

Mx  =  A  j  y  +  M|p  Y 
=' qA  a  f jCm  T  4  mT;Y  +  me  y 
^Z  =  ^  ^2  Fn,T  +  ^T,Z  +  ^E,Z 

When  staged  release  occurs,  additional  terms  required  to  define  the  full-scale  moments  and 
forces  are  given  in  Appendix  I. 
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APPENDIX  I 

DYNAMIC  EQUATIONS  OF  MOTION  MODULE 


The  flow  diagram  for  the  dynamic  equations  of  motion  module  is  shown  in  Fig.  1-1. 
Positive  directions  for  the  full-scale  forces  and  moments  are  given  in  Fig.  1-2,  and  graphic 
descriptions  of  pivoting,  rail-restricted,  and  ejector-restricted  motion  are  presented  in  Fig. 
1-3.  The  equations  for  calculating  restraining  forces  and  moments,  accelerations,  and 
velocities  are  listed  as  follows  in  blocks  as  noted  on  the  flow  diagram. 


,0  r 


Unrestricted  Motion  (Section  A) 


d  Vc 


rr\ 


Vo  5ooW 

x  vc ) 


u  =  F../m  -  qw  +  rv 


x  |  v  =  F’y  rn  — ru  +  pw 


t/* 

<0.  ■  • 

V  ' : 


,  f.  ,  .  T?  \ 

r  Any 


Fy/m  —  pv  +  qu 

- 


Uj \  vO  -t  cv  i j ^ o  r 


'-■>  r  <7  e  v  ’  -  ••••  io-.S  € 


oi-  av’bi  v/cAuy 


, 

xi 

u 

4~y-/- 

Yl 

=  [TOBQDYj 

V 

J  P 

h 

w 

-Ce\  *  \  v  & 

% 

fO  c>,  CO%-vM 

ioAy  Ccc* 

1  )--a  d. 

a  ^  f 

\./  4?  \  o  ("Y  t  ' 

IaJ  1  -S  j 

-V  y 

Ja 

dS- 

1  -x 

ou 

0^  Y 

1  iiaAvj 

r  1 1 


uj  -  ■p  i  +  v  r  X 
I  -  i  r 

-j 

c  ~  ><  L  1-  y  |  V  e-  k. 


Hv 


H, 


H- 


«x 


=  [I] 


H  x 
14  ■ 


tody-  axis)  linear 
accelerations 


Inertial  —axis 
linear  velocities 

X  X  -  7*  y 
--Tx  y  ty 

l-r  y£ 

Angular  momentum 

=  ’  ply  -  CL  1  y-/ 


■  _Z 

Ty  ?  j 

/  —  i 

X*  j 
r  Tx  g 


( +  rHy  -  qnz  • 


x  LAMMr 
M 


II Y  =  My  +  pHz  -rHx 

'  r/Uyrvyi 

lly  =  M  2  +  qll^  —  pHy 


-  pixv  1-|X7  -  Mx/£ 
-“fX'/i  +  r  x  2 

Angular  momentum 

derivatives  . 

fcss.wv-*  i A 

,x4-  X  Z-  o(~  Sy  m  nsAf'f 


J 


*-  V  a 


yy 


p 

V 

q 

-CO"1 

hy 

r 

HZ 

Body— axis  angular 
accelerations 
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Inertia  Conversion  (Section  Bl) 


•  °  s  i0 , 


XX 


1  XY 


XZ 


YY 


YZ 


zz 


XX  +m(Yo+  Zo) 


Ivv  + 

I 


‘XY 


X  Y 


I Y  7  +  m  X  Z 

a  L  o  o 


AYY 


lZZ 


Y  Z 

o  o 


(X2+Yo) 


[i] 


*xx 

-I  XY 

~  rxz 

_I  XY 

1  YY 

-I  YZ 

~I  XZ 

_I  YZ 

*ZZ 

Inertia  matrix  about 
rotation  center 


[i] 


-1 


=  inverse 


of 


[-] 


Inverse  inertia  matrix 
about  rotation  center 


Hook  Motion  from  Accelerating  Aircraft  (Section  B2) 


W.jrf's  tV'jXA"  \o-  OJJvoi-ff-t 


n 


^  J  (_1\  /  ^  'l  IV1  (■  I  ^ 


(9  =  q„  (t)  57.2958  _  ,  \ 


f  BB]  =  TME  ti/  +  (?  ,  /  ,  0,  1)  Inertial —  to —pylon  matrix 

r  ^  (glVc-inj  j  j  iClit  C-^ 

Pylon —to —  body  matrix 


[bpylon]  =  [tobody]  [bb] 

r  i  T 

[AAj  =  TME  (0,  -0p,  0,  2^  (  Y 

[ee]  =  [tobody]  [aa] 


Ot  laA  P  p»  ° J  -Zo>  1  1 

Flight— to— body  matrix 


a  Y 

O  ,  A 

0 

O 

CO 

=  [EE] 

0 

a  y 
o,  L 

a  rj 

z,p 

Resolve  A/C  (hook) 
acceleration  to  body  axis 
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Pivot  Motion  Equations  (Section  B3) 

For  all  pivot  motion,  Y0  must  be  zero;  hence,  Ixy  =  Iyz  =  0 


Rm  =  0 


V  7  =  a  7  t 

Z,p  Z,p 


vo,X 

0 

vo,Y 

=  [ee] 

0 

Vo,Z 

V7 

z,p 

Pitch  reaction  moment  always  "0" 
A/C  (hook)  velocity 


Resolve  hook  velocity  into  body  axis 


=  sin  BPYLON  (1,3)J  Store  pitch  motion  travel 

Reaction  Moments,  Pivot  Motion,  Pitch  Only  (Section  B4A) 

Set  reaction  moments  to  force  Hx  =  Hz  =  0;  hence,  p  =  r  =  p  =  r  =  0. 

R('  =  “[MX  -Zo(FY -mao,Y)  -mXo  1vo,y] 

R n  =  ~[MZ.+  XQ  (fy  -ma0,Y)  "mZo  ^o.y] 

Reaction  Moments,  Pivot  Motion,  Pitch  and  Yaw  (Section  B4B) 

Set  reaction  moment  to  force  p  =  0;  hence  p  =  0. 

Rn  =  0 

R('  =  ~^XZ'^ZZ^  [MZ  +  Xo^  Y  _mao,Y^  ~mZo  *lvo,Y 

-Mx  +  Zo[FY  -rn(ao,Y  +  rvo,X)]  +  mXo  (<tvo,Y  +  rvo,Z^ 

-  qr  (I YY  -Izz) 

Reaction  Moments,  Pivot  Motion,  Pitch,  Yaw,  and  Roll  (Section  B4C) 

R«  =  0 


R  =  0 

n 
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Rail  Motion  (Section  B5) 


XL,S 

-X 

0 

yl,s 

=  [tobody] 

-Y 

0 

’  ZL,S 

-Z 

o 

Define  hook  reference  in  inertia]  axis 


lC,H  =  XL,S+  XI  +  Rp  sin0p 


1  C,H  =  1  L,S  +  1  I 

ZC,H  =  ZL,S  +  ZI  "  Rp  ^  _  cos  V 

if  «?p  a  o),  xCH  =  xCH  +  UA(t) 


XP,H 

XC,H 

yp,h 

■  M 

yc,h 

ZP,H 

ZC,H 

X 


p,ck  =  Ap,e 


Define  hook  travel  along  pylon 


FP,X 

FX 

F  P,Y 

=  [ bpylon]' 

fy 

FP,Z 

Fz 

Resolve  full-scale  forces  into  pylon  axis 


aP,X 

a  v 

O,  A 

aP,Y 

=  [_ bpylonJ 

a  v 

O,  I 

aP,Z 

a  ry 

o,Z 

aP,X  ~ 


Resolve  hook  accelerations  into  pylon  axis 


Define  hook  accelerations 
along  pylon  X— axis 
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UB 

aP,X 

VB 

=  [bpylon] 

ap  Y 

vVB 

aP,Z 

ao,X  =  UB 
ao,Y  =  ^B 
ao,Z  =  WB 


Resolve  hook  accelerations  into  body  axis 


Hook  accelerations 


vo,  X  =  UB  J 

vq  y  =  Vg  \  Hook  velocities 

vo,Z  =  WB  \ 

Reaction  Moments,  Rail  Motion,  Translate  and  Pitch,  Side  Rail  Only  (Section  B6A) 


Z0  must  be  zero;  hence  Ixz  =  Iyz  -  0*  Reaction  moments  force,  p,  r  =  0;  hence  p 

r  =  0. 

Rm  =  0 

Rf  =  -MX-FZYoi+mYo(ao,Z-clVo,x)  +  mXoclvo>Y 

-(Txy/IyY^  [MY  “  FZX0  +  mXo  a0,Z  J" 

Rn  =  -MZ  —  FYXo  +  FXYo  +  mXo  ao, Y  _mYo  ^ao,X  +  Rvo,Z^  +  ^^XY 


Reaction  Moments,  Rail  Motion,  Translate  and  Yaw,  Bottom  Rail  Only  (Section  B6B) 

Y0  must  be  zero;  hence  Ixy  =  Iyz  =  0.  Reaction  moments  force  p,q  =  0;  hence  p  =  0 
and  q  =  0  or  constant. 

Rn  =  0 

R('  =  -Mx  +  fyZo  ~mZo  (ao,Y  +  rvo,X^  +  mXo  (rvo,Z  +  1vo,y'  +  ‘l1’  (*ZZ  _IYY^ 
-(hz^ZZ*  [MZ  +  FYXo-mXoao,Y-mZo  1vo,Y“^Xz] 

Rm  =  -MY  +  FZXo-FXZo  +  mZoK,X  -rv0;Y)  -mXo  ao,Z  -r^XZ 
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Reaction  Moments,  Rail  Motion,  Translate,  Pitch  and  Yaw  (Section  B6C) 


Either  Y0  or  Z0  must  be  zero;  hence  Iyz  =  0.  Reaction  moments  force  p  =  0;  hence 
P  =  0. 


R 

=  R 

_ 

0 

in 

n 

Rp 

=  - 

Mx  + 

FyZ 

I  o 

-f’zYo  +  ">Yo  K,Z- 

-qvo<x) 

-  (a  y-f  rv 

O  x  0,1 

o.X^ 

+ 

mXo 

(3vo,Y 

+  rvo,Z)  +  V-^ZZ  “ 

iyy) 

- 

(TXY 

'Tyy)  [ 

M  Y  —  ^  ZX  o  +  R  X  Z  o 

-mZo  (ao,X  - 

_  rvo,Y^ 

+ 

mX 

o 

S  7  — 

o ,  A 

mYo  rvo,Z  +  Rr  *XY  + 

r2^xz] 

- 

(IxZ 

;  '-I zz) 

[Mz  +  FvX„-FxY0 

_rr,Xo  ao,Y  + 

mYo  (ao,X 

+  3vo,z) 

- 

mZo 

Svo,Y  ' 

-q2TXY-qrTxz] 

Pylon  Axis  Reaction  Moments  (Section  B7) 


RP,f 

r  -i ' 

Rf 

RP,m 

=  [bpylonj 

Rm 

RP,n 

R  n 

Rail  Translation  Only  Motion  (Section  Cl) 


For  this  case,  a  rotation  center  is  undefined;  consequently,  eg  velocities  and  accelerations 
are  equal  to  hook  velocities  and  accelerations. 


R 


e  ”  “Mx  +  fyZo  “fzYo  +  mYo  K,z  -Rv0,x)  ~mZ0  K,y  +  rvo,x) 
-Mv 


R  = 

n 


-M. 


RP,f 

re 

RP.m 

=  [bpylon] 

Rm 

RP,n 

R 

n 

Pylon— axis  reaction  moments 
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P  =  o 

q  =  o 
f  =  0 


/ 

\ 


“  =  “B 
v  =  vB 

W  =  W  g 


I 

\ 


*1 

r  i ' 

u 

Yi 

=  [tobodyJ 

V 

h 

w 

Body -axis  angular  accelerations 


Body -axis  linear  accelerations 


Inertial —axis  linear  velocities 


RP,X  =  0 

Rp  Y  =  m  ao,Y-—  ^  P,Y  Pylon —axis  reaction  forces 

Rp  y  =  maoZ  —  RP,Z 


Rx 

RP,X 

ry 

=  [bpylon] 

rp,y 

Rz 

RP,Z 

Body— axis  reaction  forces 


Restrained  Motion  Rotational  Dynamics  (Section  C2) 


HX 

p 

By 

-m 

q 

Rz 

r 

Angular  momentum 
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»x  = 

=  M 

X  + 

1-..Y  - 

L  o 

Rr  + 

1 1 1 Y 

-qlB 

/  “ 

111X 

0 

(qv 

0,  Y  + 

rv 

0,  l) 

- 

niY 

(a  ry  - 

o  '  o,Z 

i 

< 

c 

4  mZ 

(a 

O  V  O 

,  Y  +  rvo 

.x> 

1*1  Y  = 

=  M 

Y  ~ 

FZXo  + 

FXZo  + 

Rn  + 

m 

P!1Z 

;  -rI1 

X  +  mX0 

(a 

0,  Z  + 

Pv 

0 ,  Y^ 

- 

mY 

o  (rvo,Z 

+  Pvo,X> 

—  m? 

'oK 

0,  X  “ 

rvo,Y> 

11  z  . 

=  M 

Z  + 

FYXo- 

FyY  + 

A  o 

«n  + 

-pH 

Y  ' 

-  mX 

0 

(a 

0, Y  ~ 

pv 

o,Z^ 

+ 

mY 

(a  y  -t 
o  v  o, A 

■  dVz.) 

-  mZ 

0  (pv 

0,  X  + 

qv 

0,  Y^ 

P 

HX 

q 

-  [T]“ 

IIY 

r 

HZ 

Body  — axis  angular  accelerations 


Restrained  Motion  Translational  Dynamics  (Section  C3) 


If  (MOTION  #3),  p  =  0 


u 

V 

w 


X1 

u 

p  i  * 

Yi 

=  [tobodyJ 

V 

k 

w 

=  a  Y  +  qZ  —  rY 

O ,  A  To  O 

=  a  v  +  —  pZ„ 

o,Y  o  r  o 


ao,Z 


Body— axis  linear  accelerations 


Inertial —axis  linear  velocities 


Reaction  Forces  (Section  C4) 


aX  =  ao,X  ~Xo  ^2  +  r2)  +  Yo  +  Z0  (4  +  Pr)  +  qvo,Z  -rvo,Y 

aY  =  ao,Y  +  Xo(f  +  Pq)  “Yo  (r2  +  P2)  +  Zo  ((Ir“P)  +  rvo,X  “Pvo,Z 
aZ  =  ao,Z  +  Xo  (Pr  _<i)  +  Yo  (P  +  V)  ~Zo  (P2  +  +  Pvo,Y  -*Ivo,X 
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Body  —axis  reactions  forces 


Pylon— axis  reaction  forces 


Body  — axis  reaction  forces 


Ejector  Plane  Restraint  Equations  (Section  C5C) 

0p  =  qp  (t)  57.2958 


[EJ  ]  =  TME  (vo  +  dp,  rj Q,  C0in  +  IR,  1) 

Inertial —to —ejector  matrix 

[ebody]  =  [tobody]  [ej] 

Ejector— to  — body  matrix 

[aa]  =  TME  (0,  -dp,  0,  2) 

[ee]  =  [tobody]  [aa] 

Flight —to  — body  matrix. 

[bB  ]  =  TME  (u0  +  dp,  Vo,  0,  1) 

Inertial —to —pylon  matrix 

[bpylon]  =  [tobody]  [bb]' 

Pylon  — to-body  matrix 

R 


max  —  Fx 


Y  Y 
ma  'y  —  P  *7 


RP,X 

RX 

rp,y 

=  [  bpylon] 

ry 

RP,Z 

RZ 

Rail  Reaction  Forces  (Section  C5B) 


R 


P,X 


RX. 

ry 

=  [bpylon] 

RP,X 

rp,y 

Rz 

RP,Z 
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i 


M 


X ,  B 


Mv  +  qr  (lyv  ~  I77) 


my,b  =  MY  +  Pr  ^ZZ  -Ixx) 
MZ,B  =  MZ  +  PI  ^XX  ~IYY^ 


Unrestrained  body  — axis  moments 


mx,e 

mx,b 

my,e 

'  -  [ebody] # 

my,b 

mz,e 

mz,b 

Unrestrained  ejector  — axis  moments 


R£,E  =  _MX,E 
Rm,E  =  0 


R 


n,E  =  _MZ,E 


Reaction  moments,  ejector  axis 


Rf 

r£,e 

R  m 

=  [ebody] 

Rm,E 

R 

n 

Rn,E 

Reaction  moments,  body  axis 


RPi 

R£ 

R  P  ,m 

="[bpylon] 

Rm 

RP,n 

R 

n 

Reaction  moments,  pylon  axis 


Hx  -  Mx  +  qr  (IYY -Izz)  +  R(/ 
Hy  =  My  +  pr  (Izz  -IXX^  +  Rm 


Angular  momentum  derivatives 


Rz  =  Mz  +  pq  (IXX  -JYY>+  Rn 
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4  -t-r 


lody  —  axis  angular  accelerations 


[ee] 


Resolve  A/C  acceleration  to  body  axis 


[eBODy]'  ao1 


Resolve  A/C  acceleration  to  ejector  axis 


fy  e  -  [  ebody]  '  fy 


Resolve  unrestrained  body— axis  forces  into 
ejector  axis 


Rx  E  =  m  aX,E  ~FX,E 
R Y.E  =  m  aY,E  “FY,E 


Reaction  forces,  ejector  axis 


ry  =  [ebody]  r. 


Reaction  forces,  body  axis 
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RP,X 

Rx 

^P,Y 

=  [ blylon]" 

ry 

Rp,z 

Rz 

Reaction  forces,  pylon  axis 


u.  =  (Fx  +  Rx)/m  —  qw  +  rv 
v  =  (F  y  +  Ryl/m  —  ru  +  pw 
w  =  (F^  +  R2;)/m  —  pv  +  qu 


Body-^axis  linear  accelerations 


x. 

u 

Yi 

=  [tobodyJ 

V 

Inertial —axis  linear  velocities 

h 

w 

Termination  Equations  (Section  C6) 


Rx  =  0 

RP,X 

o 

ii 

05 

rp,y 

Rz-  -  0 

RP,Z 

R(/  =  0 

Rp,P 

Rm  =  0 

m 

RP,m 

R  =  0 

n 

RP,n 

MOTION  =  0 


Assign  Integrator  Inputs  (Section  D) 

PI  =  u 
P2  =  v 
P3  =  w 
P4  =  Xj 
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P5  =  Yj 
P6  =  Zj 

P7  =  p 
P8  =  q 
P9  =  r 


P10  = 

TOBODY  (2,1)' 

r  — 

TOBODY  (3,1)] 

Pll  = 

'tOBODY  (3,1)" 

P  - 

'TOBODY  0,1)" 

P12  = 

"TOBODY  (1,1)' 

q  - 

'TOBODY  (2,1)" 

P13  = 

'TOBODY  (2,2)' 

r  — 

TOBODY  (3,2)" 

P14  = 

TOBODY  (3,2)' 

P  - 

'TOBODY  (1,2)’ 

P15  - 

TOBODY  (1,2)' 

q  - 

'TOBODY  (2,2) 

P16  = 

TOBODY  (2;3)" 

r  — 

TOBODY  (3,3)  _ 

P17  = 

TOBODY  (3,3)~ 

P  - 

'TOBODY  (1,3)" 

P18  = 

TOBODY  (1,3)' 

q- 

"TOBODY  (2,3)] 

P19  =  uB 
P20  =  v  B 
P21  =  vvB 
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MOTION  OPTIONS 

0  UNRESTRAINED 

1  PIVOT;  PITCH  ONLY 

2  PIVOT;  PITCH  AND  YAW 

3  PIVOT;  PITCH,  YAW  AND  ROLL 

4  RAIL;  TRANSLATE  ONLY 

5  RAIL;  TRANSLATE  AND  PITCH 

6  RAIL;  TRANSLATE  AND  YAW 

7  rail;  TRANSLATE,  PITCH  AND  YAW 

0  TRANSLATE, ROTATE  IN  EJECTOR  PLANE 


Figure  1-1.  Flow  diagram  of  the  dynamics  module. 
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TOP  VIEW 


Figure  1-2.  Positive  directions  of  full-scale  forces  and  moments. 
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CONTROL  PARAMETER  IS  MOTION 

INITIAL  RELEASE  FINAL  RELEASE 

MOTION  TYPE  MOTION  CRITERION  CRITERION 


0 

UNRESTRICTED 

— 

— 

1 

PIVOT;  PITCH  ONLY 

A0r 

Rp,Z 

2 

PIVOT;  PITCH  AND  YAW 

A0r 

Rp,e 

3 

PIVOT;  PITCH,  YAW,  ROLL 

A0r 

Rp.z 

4 

RAIL;  TRANSLATE  ONLY 

xp,i 

xp,i 

5 

RAIL;  TRANSLATE  AND  PITCH 

XP,I 

XP,2 

6 

RAIL;  TRANSLATE  AND  YAW 

Xp,l 

XP,2 

7 

RAIL;  TRANSLATE,  PITCH  AND  YAW 

XP,I 

XP,2 

8 

TRANSLATE,  ROTATE  IN  EJFCTOR  PLANE 

PIVOT  MOTION  {OPTIONS  1,2,3) 

RESTRICTION:  Yo=0 

EJECT 

EJECT 

INITIAL  IF  (Afl  -A0r  )  <0,  CHECK  RF>2 

FINAL  IF  Rpz  <  0,00  TO  UNRESTRICTED  MOTION 


a.  Pivot  motion  (options  1  through  3) 

Figure  1-3.  Graphic  descriptions  of  staged  release  options. 
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RESTRICTIONS:  4 

5 

6 
7 


NONE 

SIDE  RAIL  ONLY (Zo  =0) 

BOTTOM  RAIL  ONLY  (Yo  =  O) 

SIDE  OR  BOTTOM  RAIL(Zo  OR  Yo  MUSTS 0) 


NOTES:  o)  OPTIONS  4  —  7,  TRANSLATE  ONLY  FOR  AFT  HOOK  TRAVEL  LESSTHAN  Xpj 
b)  OPTIONS  5—7,  ANGULAR  MOTION  (AS  DESCRIBED)  IN  ADDITION  TO 

TRANSLATION  FOR  AFT  HOOK  TRAVEL  GREATER  THAN 
Xpj  BUT  LESS  THAN  Xp<2 

C)  UNRESTRICTED  MOTION  FOR  AFT  HOOK  TRAVEL  GREATER  THAN 
Xpj  (OPTION  4),  GREATER  THAN  XP)2  (OPTIONS  5—7  ) 


HOOK  FWD 


+CN 


Yo  . 


CENTER  OF 
ROTATION 


+  2B 

REAR  VIEW 

XR=  0 


■+yb 


+yb 


+cN 


CENTER  OF 
ROTATION 


J 


■  Ho 


+  €f 


REAR  VIEW 
Ir - -90  DEG 


b.  Rail  motion  (options  4  through  7) 
Figure  1-3.  Continued. 
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TRANSLATE,  ROTATE  ONLY  IN  PLANE  OF  EJECTORS  (OPTION  6  ) 

ASSUMPTION  •  MOTION  ABOUT  eg  ,  NO  INERTIA  TRANSFER  REQ'D 
RESTRICTION:  IXY  =  Ixz  =  IY2  =  0 


STORE  IS  RESTRAINED  TO  TRANSLATION  AND  ROTATION 
IN  THE  PLANE  OF  THE  EJECTORS  DURING  EJECTOR  ACTION 


RELEASE  CRITERIA:  IF  EJECT  =0,  GO  TO  UNRESTRAIN  ED  MOTION 


c.  Ejector  plane  motion 
Figure  1-3.  Concluded. 
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APPENDIX  J 

OUTPUT  PROCESSING  EQUATIONS 

Calculate  Store  Coordinates  With  Respect  to  the  Flight-Axis  System  Origin  (See  Fig.  10). 

ep  =  qp  (t)  57.2958 

For  unaccelerated  flight  ( 0P  =  0) 

xt  -  X, 


For  accelerated  flight  (0P  =£  0) 
[AA]  =  TME  (0,  -6  ,  0,  2) 

.  XC  =  XI  +  RPsin0p  +  UAt 

Yc  =  Yj 

ZC  =  Zl  -Rp  (]  ~cos0p) 


X 

t 

i  ' 

XC 

Yt 

l 

< 

< 
- 1 

II 

Yc 

Z 

t 

zc 

y  ' 


Calculate  Store  Attitude  With  Respect  to  the  Free-stream  Wind  Vector  (See  Fig.  15). 

tan-1  [Zj.'tX,  +  UA)] 


Th  -  m  -tan  1  {^l/[(XI  +  UA> 

(0^  =  CO  J 


cos  v\  —  Z T  sin  v 


>.]} 
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APPENDIX  K 

CTS  CLOSED-LOOP  POSITIONING  MODULE  EQUATIONS 

The  flow  diagram  for  the  CTS  closed-loop  positioning  module  is  shown  in  Fig.  K-l,  and 
the  equations  used  for  calculating  CTS  rig  coordinates  and  attitudes  follow. 

CTS  Rig  Angles*  (Section  A) 


> 

> 

I  =  TME  (i^,  T7t,  o,  1) 

Store  attitude  with  respect  to 

J  1  '  1 

free— stream  wind  vector 

[bb" 

j  =  TME  (Au,  A r],  0,  l) 

[dd 

]  -  [bb]'[aa] 

Account  for  deflections 

[bb 

]  =  TME  (<AS,  ds,  0,  2) 

[cc 

1  ■  [bb]'[dd] 

Account  for  sting  bend  angles 

7r 

=  sin-1  [cc  (1,2)] 

1 

l'R 

=  sin  —  1  [  — CC  (1,3) ,/cos  ] 

CTS  rig  angles  in  pitch,  yaw,  roll  sequence 

CUR 

ll 

£ 

l 

£> 

e 

\ 

CTS  Rig  Linear  Coordinates*  (Section  B) 

[AA]  =  TME  ,(,  cu.,  1) 


axr,i 

0 

ayr.i 

i"<~1 

< 

1 - 1 

II 

Ym 

aZr,i 

Z,„ 

*See  Fig.  15. 

*See  Figs.  A-l  and  A-2. 
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[aa]  =  TME  (./R,  ,R,  0,  1) 


bb]  = 

TME  (^s,  0S,  0,  2) 

cc]  = 

[bb]  [aa 

AXR,2 

0 

AYR,2 

-  [ccT 

0 

AZR,2 

~dR 

AXR,3 

^1,R 

"AYR,3 

-  [**]' 

0 

AZR,3 

0 

[BB]  =  TME  (^s,  0,  0,  2) 

[cc]  =  [bb][aa] 


AXR,4 

0 

AYR,4 

-  [cc] 

^2,R 

AZR,4 

0 

[aa]  =  TME  (vR,  0,  0,  1) 


AX  R ,  5 

1-  ”1' 

3 

AYR,5 

< 

II 

0 

AZR,5 

0 
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AXr 

AYr 

XZr 


5 

V  AX 

n=  1 


R,n 


2  AY 

n  =  1 


R,n 


l  AZp 

n  =  1  R>n 


Distance  from  CTS  pitch  center  to  store  eg 


Xt,R  =  Xtd2X) 
Yt,R  =  Y,  (12A) 

Zt,R  =  Zt  <12X) 


Convert  from  ft,  full  scale,  to  in.,  model  scale 


Xr  -  XTP  o  _AXR -AXL  +  xt  R  ' 

YR  =  VTP,o  ~AYR ~AYL  +  Yt  R  +  AYC  V  CTS  rig  coordinate 
Z  R  =  ZTP,o  ~  AZR  ~AZL  +  Zt,R  +  AZC  1 
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READ  TOUCH  FILE 


|  CLEAR  ERROR  FLAG  | 

EXIT 


i - J  •: 

NO 

i  SECT  A 

i 

CALCULATE  NEW 

RIG  ANGLES 

'  SECT  B 

CALCULATE  NEW 

RIG  POSITIONS 

_NO_ 


CONVERT  RIG  POSITIONS 
AND  ATTITUDES  TO 
9AC  COUNTS 


CALCULATE  DIFFERENCES 
IN-  COMMAND  AND  SET 
POSITIONS  AND  ANGLES 


ERROR  FLAG  =  TRAVEL  LIMIT  IN  FLOW  DIAGRAM 
OF  TRAJECTORY  GENERATION  CLOSED- LOOP  PROGRAM 


Figure  K-1.  Flow  diagram  of  the  CTS  closed-loop 
positioning  module. 
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APPENDIX  L 

NONROLLING  STING  APPLICATIONS 


For  trajectory  tests  which  use  a  nonrolling  sting  support,  roll  commands  which  are 
generated  in  the  simulation  obviously  cannot  be  executed  by  the  CTS  rig.  Therefore,  if  not 
restrained  in  some  manner,  the  simulated  roll  position  (as  determined  from  trajectory 
calculations)  and  actual  model  roll  position  could  possibly  digress  to  a  point  where 
trajectory  results  are  not  representative.  For  the  Tunnel  4T  trajectory  program,  the  restraint 
chosen  was  to  set  the  simulated  roll  position  equal  to  the  actual  model  roll  position  for  the 
pass  through  the  conversion  module  of  the  trajectory  equations  made  just  after  input 
processing.  For  stores  with  similar  planforms  in  the  Xb-Zb  and  Xb-Yb  planes,  application  of 
this  restriction  would  probably  have  only  minor  effects  on  the  trajectory  outcome  since  the 
roll  attitude  would  not  be  expected  to  significantly  affect  aerodynamic  loading.  For  stores 
with  dissimilar  planforms  in  the  Xb-Zb  and  Xb-Yb  planes,  effects  of  this  restriction  on 
trajectory  results  are  potentially  more  significant  since  aerodynamic  loading  would  be 
expected  to  change  with  roll  attitude.  Equations  used  to  apply  this  restraint  are  listed  as 
follows: 


NONROLL 


Options 


0 

1 

2 

3 


Roll  capability 

No  roil  capability,  set  simulated  roll  position  equal  to 
actual  model  roll  position  each  pass  through  the 
conversion  module  after  input  processing 

Same  as  option  1,  but  use  instead  for  0-  or  6 -in.  offset 
roll  mechanisms  with  no  roll  capacity  to  impose  additional 
sting  moment  limitations 

Same  as  option  1  except  set  Cfir=  0  in  trajectory 
calculations 


Calculate  Roll  Position  As  Set 

M 

[bb] 

[dd] 

[bb] 

[aa] 

6>TR 


TME  (iar,  77R,  0,  1) 


TME  (0S,  0gt  0,  2) 

[bb]  [aa] 

TME  (A„,  A rj,  0,  1) 

[bb]  [dd] 

tan  ~  1  [-AA(3,2)/AA(2,2)J 
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Calculate  Difference  Between  Actual  and  Simulated  Roll  (Euler  Sequence) 

[no]  =  [tobody]  [aa] 

A^6tr  =  tan-1  [rO(2,3)/RO(3,3)  ] 

Calculate  Roll  Difference  (Pitch,  Yaw,  Roll  Sequence) 

A&JTR  =  I  wl  ~  WTR  I 

If  any  of  following  checks  is  true,  bypass  the  remainder  of  the  calculations. 

If  (Aojtr  <  0.25) 

If  (NPASS  <  5) 

If  (ROLFLG  ='  1) 1 

If  not,  continue: 

Set  Simulated  Roll  Equal  to  Actual  Model  Roll 

“I  ~  WTR 

Update  Inertial-to-Body  Matrix 

■  -  [tobody]  =  TME  (i/,,  1J,,  <0,,  1) 

<f>{  =  tan-1  [TOBODY  (2,3). 'TOBODY  (3,3)] 

Update  Integrator  Outputs 

Pll  =  TOBODY  (2,1) 

P12  =  TOBODY  (3,1) 

P14  =  TOBODY  (2,2) 

P15  =  TOBODY  (3,2) 

P17  =  TOBODY  (2,3) 

PI  8  =  TOBODY  (3,3) 

*The  flag  (ROLFLG)  is  set  equal  to  one  after  the  pass  through  the  trajectory  calculations  in  input  processing 
and  is  set  equal  to  zero  just  after  the  pass  counter  is  updated  in  output  processing.  This  procedure  prevents  the 
integrator  outputs  from  being  manipulated  when  the  trajectory  modules  are  called  during  the  integration  process. 
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APPENDIX  M 
MATRIX  DEFINITIONS 

Notation  of  Terms 


a(l,l) 

a(l ,2) 

a(l,3) 

a(2,l) 

a(2,2) 

a(2,3) 

a(3,l) 

a(3,2) 

a(3,3) 

For  Pitch,  Yaw,  Roll  Sequence 


TME  (a,  Op  a9,  1)  = 


cos  a  cos  a , 


sin  a 2  sm  a 


-cos  <*2  cos  ol  sin  a  j 


sina2  cos  a  sina^ 


cos 


sin  a 


sin  an 


cos  C2  cos  a  j 


—sin  a2  cos  a  j 


—sin  a  cos  a 


.1 


cosa2  sina  smaj 
+  cos  a  sin  a0 


cos  a 2  cos  a 
—sin  a 2  sin  a  j  sin  a 


For  Euler  Sequence 


TME  (a, 


a,,  a 


,2)  = 


cos  a  cos  a 


sin  (*2  sin  a  i  cos  a 
—cos  a9  sin  a 


cos  «2  sin  a  j  cos  a 
+  sin  a2  sin  a 


cos  a  ^  sin  a 


sin  a 2  sin  aj  sin  a 
+  cos  a9  cos  a 


cos  a2  sin  a  2  sin  a 
—  sin  a 2  cos  a 


—sin  a1 


sin  cos  a 2 


cos  a9  cos  a  j 
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A/C 

an>hn 


ao,X>ao,Y> 

ao,Z 


aX>aY>aZ 


az,p 

BPYLON 


Ca.oiCn 

Cy.o 


,oj 


(Ca,o)o,(CN  ,o)oi 
(Cy,o)o 


(Ca,o) 


max 


Qa.uCn.Cy 


C-A.t.T  >Cn,T  > 

Cy,t 


^'A,t,X!  Cn,x. 
Cy,x 


Cjdf>Cjdm, 

Cjdn 


NOMENCLATURE 

Store  full-scale  and  model-scale  reference  areas,  respectively,  ft2 
Aircraft  model  designation 

Fifth-order  polynomial  curve  fit  coefficient  values  which  define 
the  thrust  force  as  a  function  of  time  (see  Appendix  F) 

Accelerations  of  the  store  rotation  center  (if  not  coincident  with 
the  mass  center),  positive  in  the  positive  Xb,  Yb,  and  Zb  directions, 
respectively,  ft/sec2 

Accelerations  of  the  store  mass  center,  positive  in  the  positive  Xb, 
Yb,  and  ZB  directions,  respectively,  ft/sec2 

Acceleration  of  the  aircraft  in  the  Zp  direction,  ft/sec2  (see  Fig.  10) 

Notation  for  the  matrix  (3x3)  which  converts  pylon-axis  quantities 
to  body-axis  quantities 

Current  values  of  the  external  input  axial-force,  normal-force,  and 
side-force  coefficients,  respectively 

Initial  values  of  the  external  input  axial-force,  normal-force,  and 
side-force  coefficients,  respectively 

Maximum  value  of  the  external  input  axial-force  coefficient  for 
ramp  axial-force  coefficient  calculations  (see  Appendix  D) 

Store  measured  axial-force,  normal-force,  and  side-force 
coefficients,  positive  in  the  negative  Xb,  negative  Zb,  and  positive 
Yb  directions,  respectively 

Sum  of  the  aerodynamic  force  coefficient  contributions  acting  on 
the  full-scale  store,  positive  in  the  negative  Xb,  negative  Zb,  and 
positive  Yb  directions,  respectively  (see  Appendix  E) 

Extrapolated  (or  measured)  values  of  the  axial-force,  normal- 
force,  and  side-force  coefficients,  respectively 

Jet  roll-damping,  pitch-damping,  and  yaw-damping  coefficients, 
respectively,  ft-sec 
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Cf,Cm,Cn 


C^oiCnijOJ 

Qi.o 

(Cf,o)oj(Cm,0)0> 

(Cn,o)o 

cfp,cmq,cnr 


Cf.T.C  m,  T,Cn,T 


Cf,xiCmiX5 

Cn,x 

ACn 


COEF,COEFI 

CONFIG 
CON  SET 


CiB,C2B 


Cn>dn 


DATE 


Store  measured  rolling-moment,  pitching-moment,  and  yawing- 
moment  coefficients,  respectively;  the  positive  vectors  are 
coincident  with  the  positive  Xb,  Yb,  and  Zb  directions 

Current  values  of  the  external  input  rolling-moment,  pitching- 
moment,  and  yawing-moment  coefficients,  respectively 

Initial  values  of  the  external  input  rolling-moment,  pitching- 
moment,  and  yawing-moment  coefficients,  respectively 

Store  roll-damping,  pitch-damping,  and  yaw-damping  derivatives, 
respectively,  per  radian 

Sum  of  the  aerodynamic  moment  coefficient  contributions  acting 
on  the  full-scale  store;  the  positive  vectors  are  coincident  with  the 
positive  Xb,  Yb,  and  Zb  axes,  respectively  (see  Appendix  E) 

Extrapolated  (or  measured)  values  of  the  rolling-moment, 
pitching-moment,  and  yawing-moment  coefficients,  respectively 

Difference  in  successively  measured  values  of  normal-force 
coefficient 

Current  and  initial  values  of  the  offset  coefficient  module  control 
flag 

Aircraft  model  configuration  designation 

Run/point  number  of  constants  set  used  in  data  reduction 

Difference  in  pressure  coefficient  between  probe  orifices  1  and  3, 
positive  for  positive  e,  (pSji  -  ps,3)/qL 

Difference  in  pressure  coefficient  between  probe  orifices  2  and  4, 
positive  for  positive  a,  (pS)4  -  ps,2)/qL 

Break  points  for  forward  and  aft  ejector  force  polynomial  curve 
fits,  respectively,  sec  or  ft  (see  Appendix  G) 

Fifth-order  polynomial  curve  fit  coefficient  values  which  define 
the  forward  ejector  force  as  a  function  of  time  or  stroke,  sec  or  ft 
(see  Appendix  G) 

Calendar  time  at  which  data  were  recorded 


115 


AE  DC-TR-79-1 


Slope  of  the  external  input  axial-force  coefficient  for  ramp  axial- 
force  coefficient  calculations  (see  Appendix  D) 

Vertical  offset  of  the  CTS  support  sting,  in.  (see  Fig.  A-2) 


dCA,0/dt 


dR 


EE 

EJECT, EJECTI 
E 1 FL  AG  ,E2FLAG 
^n>fn 


FA,g>FN,g» 

Fy.g 

Fa,i*Fn»Fy 


Fe.XsF^y, 

Fe,z 


Notation  for  the  matrix  (3x3)  which  converts  flight-axis  quantities 
to  body-axis  quantities 

Current  and  initial  values  of  the  ejector  module  control  flag 

Forward  and  aft  ejector  cutoff  control  flags,  respectively 

Fifth-order  polynomial  curve  fit  coefficient  values  which  define 
the  aft  ejector  force  as  a  function  of  time  or  stroke,  sec  or  ft  (see 
Appendix  G) 

Total  forces  measured  by  the  store  balance,  positive  in  the  negative 
Xr,  negative  ZB,  and  positive  YB  directions,  respectively,  lb 

Measured  aerodynamic  forces  acting  on  the  store  model,  positive 
in  the  negative  XB,  negative  ZB,  and  positive  YB  directions, 
respectively, lb 

Components  of  the  ejector  force  acting  on  the  store,  positive  in  the 
positive  XB,  YBi  and  ZB  directions,  respectively,  lb 


Fei.Fk 


Forward  and  aft  ejector  forces,  respectively,  lb 


Ft,X)Ft,Y) 

Ft,z 

Fx»Fy,Fz 


f 

g 


Components  of  the  thrust  force  acting  on  the  store,  positive  in  the 
positive  XB,  Yb,  and  ZB  directions,  respectively,  lb 

Components  of  total  force  acting  on  a  free-falling  store,  positive  in 
the  positive  XB,  YB,  and  ZB  directions,  respectively,  lb 

Notation  for  derivatives  stored  in  the  derivative  history  file  (see 
Fig.  11) 

Acceleration  of  gravity,  ft/sec2 


Hx,Hy,Hz  Components  of  the  angular  momentum;  the  positive  vectors  are 

coincident  with  the  positive  XB,  YB,  and  ZB  axes,  respectively,  ft- 
lb-sec 

h  Simulated  pressure  altitude,  ft 

I  Notation  for  the  inertia  matrix  of  a  store  which  is  calculated  about 

the  mass  center  (see  Appendix  I) 
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Ir 

ITR 
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^ints^int 

KPASS 

KY,Y,KY,n 


Notation  for  the  inverse  of  the  inertia  matrix  (I) 

Notation  for  the  inertia  matrix  of  a  store  which  is  calculated  about 
a  point  (rotation  center)  other  than  the  mass  center  (see 
Appendix  I) 

Notation  for  the  inverse  of  the  inertia  matrix(f) 

Pitch  and  yaw  incidence  angles  of  the  store  longitudinal  axis  at 
carriage  with  respect  to  the  aircraft  longitudinal  axis,  positive  nose 
up  and  nose  to  the  right,  respectively,  as  seen  by  the  pilot,  deg 

Roll  incidence  of  the  store  ZB  axis  at  carriage  with  respect  to  the 
aircraft  plane  of  symmetry,  positive  for  clockwise  roll  looking 
upstream,  deg 

Number  of  iterations  in  the  flow-field  calculations  required  for  the 
difference  in  successive  calculations  of  0T  to  be  less  than  0. 1  deg 

Full-scale  moments  of  inertia  about  the  store  Xb,  Yb,  and  Zb  axes, 
respectively,  slug-ft2,  and  referenced  to  the  store  mass  center 

Full-scale  moments  of  inertia  about  the  store  Xb,  Yb,  and  Zb 
directions,  respectively,  slug-ft2,  but  referenced  to  a  point 
(rotation  center)  other  than  the  mass  center 

Full-scale  products  of  inertia  in  the  store  Xb-Yb,  Xb-Zb,  and  Yb- 
ZB  planes,  respectively,  slug-ft2,  and  referenced  to  the  store  mass 
center 

Full-scale  products  of  inertia  in  the  store  Xb-Yb,  Xb-Zb,  and  YB- 
ZB  planes,  respectively,  slug-ft2,  but  referenced  to  a  point  (rotation 
center)  other  than  the  mass  center 

Current  and  initial  values  of  the  parameter  which  determines  the 
number  of  integrator  passes  per  data  cycle 

Multiple-pass  integration  counter 

Linear  deflection  constants  of  the  store  model  support  sting  in  the 
Yb  direction  and  measured  at  the  balance  center,  in. /lb  and 
in. /in. -lb,  respectively 
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Linear  deflection  constants  of  the  store  model  support  sting  in  the 
Zb  direction  and  measured  at  the  balance  center,  in. /lb  and 
in. /in. -lb,  respectively 

Angular  deflection  constants  of  the  store  model  support  sting  in 
the  XB-ZB  plane,  deg/lb  and  deg/in. -lb,  respectively 

Angular  deflection  constant  of  the  store  model  support  sting  in  the 
Yb-Zb  plane,  deg/in. -lb 

Angular  deflection  constants  of  the  store  model  support  sting  in 
the  Xb-Yb  plane,  deg/lb  and  deg/in. -lb,  respectively 

Store  model  scaling  change  factor  (normally  equals  1) 

Wing' identification  print  control  constants 

Store  full-scale  reference  dimensions  for  pitching-moment, 
yawing-moment,  and  rolling-moment  coefficients,  respectively,  ft 

Store  model-scale  reference  dimensions  for  pitching-moment, 
yawing-moment,' and  rolling-moment  coefficients,  respectively,  in. 

Physical  dimensions  of  the  CTS  support  sting  used  in  model 
positioning  .and  clearance  calculations,  in.  (see  Fig.  A-2) 


Mb 


Nominal  test  Mach  number 


Me,x»Me,y.  Components  of  ejector  moment;  the  positive  vectors  are 

Mb,z  coincident  with  the  positive  Xb,  Yb,  and  ZB  axes,  respectively, 

ft-lb 

Ml  Local  Mach  number  calculated  from  probe  ps  measurements 

Mf,Mm,Mn  Measured  aerodynamic  moments  acting  on  the  store  model;  the 

positive  vectors  are  coincident  with  the  positive  XB,  Yb,  and  ZB 
directions,  respectively,  in. -lb 

Mfig,Mmig,  Total  moments  measured  by  the  store  balance;  the  positive  vectors 

Mn>g  are  coincident  with  the  positive  XB,  YB,  and  ZB  axes,  respectively, 

in. -lb 

MODE  Parameter  which  defines  type  of  CTS  operation 


MOTION,  HOLDI  Current  and  initial  values  of  staged  separation  control  parameter 
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Mt,Z 
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M 

00 

m 

NDX 

NOROLL 

NPASS 

NX 

NXP 

Nz 

P,P1  -  P50 

P,P1  -  P50 

PASS 

PN 

POST 

P.q.r 


Pa 

Po>Qo>ro 


Ps,l  ~ *  Ps,4> 

Pp,5 


Components  of  the  thrust  moment  acting  on  the  store;  the  positive 
vectors  are  coincident  with  the  positive  Xb,  Yb,  and  Zb  axes, 
respectively,  ft-lb 

Components  of  the  total  moment  acting  on  a  free-falling  body;  the 
positive  vectors  are  coincident  with  the  positive  XB,  YB,  and  Zb 
axes,  respectively,  ft-lb 

Wind  tunnel  free-stream  Mach  number 

Store  mass,  slugs 

Sequential  indexing  number  for  referencing  data  obtained  during  a 
grid  set;  indexes  for  each  position  in  the  set 

CTS  rig  roll  control  parameter 

Integrator  pass  counter 

Number  of  extrapolators 

Extrapolator  pass  counter 

Aircraft  “g”-loading  factor 

Integrator  output  designation 

Integrator  input  designation 

Data  cycle  pass  counter 

Data  point  number 

Launch/postlaunch  control  parameter 

Store  angular  velocity  about  the  Xb,  Yb,  and  Zb  axes,  respectively; 
the  positive  vectors  are  coincident  with  the  positive  XB,  Yb,  and  Zb 
axes,  rad/sec 

Static  pressure  at  the  simulated  altitude,  psfa 

Initial  values  of  the  store  angular  velocity  about  the  XB,  YB,  and 
Zb  axes,  respectively,  rad/sec 

Measured  pressures  for  probe  orifices  1  through  5,  respectively, 
psfa 
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Pt 

Pt,p 


P 


ao 


P5 


Ps.p 

Qa 

Ql 

qP 


Recc 


R  i  ;  R  Ili  ■/  R  ,1 


ROLFLG 

Rp,f>Rp,m>  Rp,n 


Rp,X>Rp,Y>Rp,Z 


Rp 

RUN 


Rx.Ry.Rz 


fiRi,«Rs 


Wind  tunnel  total  pressure,  psfa 

Probe  measured  free-stream  total  pressure  corrected  for  local 
Mach  number,  psfa 

Wind  tunnel  free-stream  static  pressure,  psfa 

Ratio  of  the  average  of  the  four  static  pressures  and  the  probe  total 
pressure,  (ps,i  +  ps>2  +  ps,3  +  Ps,4)/4(pp,5) 

P5  ratio  corrected  for  probe  attitude,  equivalent  value  for  0T  =  0 

Dynamic  pressure  at  the  simulated  altitude,  psf 

Local  dynamic  pressure,  psf  (from  probe  measurements) 

Pitch  rate  of  the  aircraft  during  a  pullup/pushover  maneuver, 
rad/sec  (see  Fig.  10) 

Wind  tunnel  free-stream  dynamic  pressure,  psf 

Wind  tunnel  free-stream  unit  Reynolds  number,  millions. per  ft 

Full-scale  body-axis  restraining  moments  about  the  pivot  (rotation 
center);  the  positive  vectors  are  coincident  with  the  positive  XB, 
Yb,  and  ZB  directions,  respectively,  ft-lb 

Roll  simulation  control  parameter  (see  Appendix  L) 

Full-scale  pylon-axis  pivot  (rotation  center)  restraining  moments; 
the  positive  vectors  are  coincident  with  the  positive  Xp,  YP,  and 
ZP  directions,  respectively,  ft-lb 

Full-scale  pylon-axis  pivot  (rotation  center)  restraining  forces, 
positive  in  the  positive  Xp,  Yp,  and  Zp  directions,  lb 

Effective  rotation  arm  of  the  aircraft  during  a  pullup/pushover 
maneuver,  ft  (see  Fig.  10) 

Data  set  identification  number 

Full-scale  body-axis  pivot  (rotation  center)  restraining  forces, 
positive  in  the  positive  XB,  YB,  and  ZB  directions,  lb 

Differential  balance  readings  because  of  model  weight  at  wind-off 
conditions  for  balance  gages  1  and  5,  counts 
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Wind  tunnel  specific  humidity,  lbm  H2O  per  lbm  air 

Current  and  initial  values  of  the  data  cycle/integration  interval 
increase  control  parameter 

Store  model  designation 

Configuration  indexing  number  used  to  correlate  data  with  the  test 
log;  survey  may  be  used  to  identify  all  or  portion  of  a  grid  set 

Static  temperature  at  the  simulated  altitude,  °R 

Wind  tunnel  dewpoint  temperature,  °R 

Alphanumeric  notation  for  referencing  a  specific  test  program  in  a 
specific  test  unit 

Current  and  initial  values  of  the  thrust  module  control  flag 

Time  at  which  data  were  recorded  (hr/min/see) 

Notation  for  current  and  initial  values  of  the  matrix  (3x3)  which 
converts  inertial-axis  quantities  to  body-axis  quantities 

Number  which  identifies  a  particular  set  of  trajectory 
constants/test  conditions 

Wind  tunnel  total  temperature,  °R 

Wind  tunnel  free-stream  static  temperature,  °R 

Trajectory  time  from  the  instant  of  store  release  from  the  aircraft, 
sec 

Data  acquisition  time  increment,  sec 
Trajectory  integration  time  increment,  sec 

Time  delay  before  initiation  of  the  external  input  ramp  axial  force, 
sec  (see  Appendix  D) 

Time  delay  before  initiation  of  the  thrust  force,  sec  (see 
Appendix  F) 

Internally  calculated  time  delay  parameter,  sec 
Initial  value  of  trajectory  time,  sec 
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tt 

tt,Cldt,C2 

t2 


UA 

Ur 


u,v,w 


Ub,VB,WB 


u0,v0,w0 


Vl 

Vx,VY,Vz 


VxY.Vxz,  VYZ 


V0)XiV0,Y. 

V0,Z 


WA,WN,WY 


WING 


Time  from  thrust  initiation,  sec 

Thrust  module  event  designation  parameters,  sec  (see  Appendix  F) 

Break  point  for  thrust  force  polynomial  curve  fits,  sec  (see 
Appendix  F) 

Velocity  of  the  aircraft  at  the  simulated  altitude,  ft/sec 

Total  velocity  of  the  full-scale  store  with  respect  to  a  space-fixed 
point,  ft/sec 

Velocities  of  the  full-scale  store  relative  to  the  origin  of  the  flight- 
axis  system,  positive  in  the  positive  XB,  YB,  and  ZB  directions, 
respectively,  ft/sec 

Velocities  of  the  hook  (rotation  center)  relative  to  the  origin  of  the 
flight-axis  system,  positive  in  the  positive  XB,  YB,  and  ZB 
directions,  respectively,  ft/sec 

Initial  values  of  the  velocity  of  the- full-scale  store  relative  to  the 
origin  of  the  flight-axis  system,  positive  in  the  positive  XB,  YB,  and 
ZB  directions,  respectively,  ft/sec 

Local  velocity,  ft/sec  (from  probe  measurements) 

Velocity  components  relative  to  a  space-fixed  axis  and  parallel  to 
the  XB,  Yb,  and  ZB  axes,  positive  in  the  -XB,  YB,  and  -ZB 
directions,  respectively,  ft/sec 

Velocity  components  relative  to  a  spaced-fixed  axis  in  the  body- 
axis  XB-YB,  XB-ZB,  and  YB-ZB  planes,  respectively,  ft/sec 

Velocity  of  the  aircraft  in  the  Zp  direction,  ft/sec 

Wind  tunnel  free-stream  velocity,  ft/sec 

Components  of  rotation  center  velocity,  positive  in  the  positive 
XB,  YB,  and  ZB  directions,  respectively,  ft/sec 

Model  weight  tares  along  the  XB,  ZB,  and  YB  axes,  respectively,  lb 

Location  of  store  launch  position 

Store  full-scale  weight,  lb 
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Wx,Wy,Wz 


WX,WY,WZ 


X,Y,Z 


AXae 

Xbf 

Xc.Yc.Zc 


Xcg 

Xei»XE2  . 


Xfe 


Xi,YIfZ, 


X^oj  Yii0, 

Zl,o 


axl,ayl,azl 


X 
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Components  of  the  full-scale  store  weight,  positive  in  the  positive 
Xj,  Yj,  and  Zi  directions,  respectively,  lb 

Components  of  the  full-scale  store  weight,  positive  in  the  positive 
Xb,  Yb,  and  Zb  directions,  respectively,  lb  (including  effects  of 
simulated  dive  or  bank  angles) 

Separation  distance  of  the  store  eg  from  the  flight-axis  system 
origin  in  the  positive  XF,  YF,  and  ZF  directions,  respectively,  ft, 
full  scale 

Distance  between  forward  and  aft  ejector  pistons,  ft,  full  scale 

Distance  from  the  store  model  nose  to  the  balance  face,  in. 

Separation  distance  of  the  store  eg  from  the  flight-axis  system 
origin  in  the  positive  Xi,  Yi,  and  Zj  directions,  respectively,  ft,  full 
scale 

Axial  distance  from  the  store  nose  to  the  eg  location,  ft,  full  scale 

Full-scale  distances  from  store  center  of  gravity  to  line  of  action  of 
•forward  and  aft- ejector  forces  measured  along  XB  axis,  positive  if 
force  is  forward  of  center  of  gravity,  ft,  full  scale 

Axial  distance  from  the  store  nose  to  the  forward  ejector  piston, 
ft,  full  scale 

Separation  distance  of  the  store  eg  from  the  inertial-axis  system 
origin  in  the  positive  Xj,  Y],  and  Z\,  directions,  respectively,  ft, 
full  scale 

Positions  of  the  store  eg  with  respect  to  the  carriage  position  at 
trajectory  initiation,  positive  in  the  positive  Xj,  Y],  and  Zi 
directions,  respectively,  ft,  full  scale 

Linear  deflections  of  the  CTS  model  support  sting,  positive  in  the 
positive  Xj,  Yi,  and  Zi  directions,  respectively,  in. 

Axial  distance  from  the  store  nose  to  the  eg  location,  in.,  model 
scale 


AXm>Cg,AXniCg  Axial  distances  from  the  store  center  of  gravity  to  the  pitching- 

moment  and  yawing-moment  reference  centers,  respectively, 
positive  in  the  positive  Xb  direction,  ft,  full  scale 
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X0,Y0,Z0 


Xp,Yp,Zp 


Xp,CK 


XpiTC.YpiTC. 

ZpiTC 

Xpt, 

Xp,2 


Xr,Yr,Zr. 


AXr,AYr, 

AZr 


Xref>Yref, 

Zref 


Xt,Y„Zt 


Axial  distances  from  the  balance  face  to  the  balance  electrical 
center  in  the  pitch  and  yaw  planes,  respectively,  in. 

Axial  distances  from  the  balance  electrical  center  to  the  store 
pitching-moment  and  yawing-moment  reference  centers, 
respectively,  positive  in  the  positive  XB  direction,  in.,  model  scale 

Distances  from  the  pivot  (rotation  center)  to  the  store  center  of 
gravity  along  the  XB,  Yb,  and  Zr  axes,  respectively,  positive  in  the 
positive  Xr,  Yb,  and  ZB  directions,  ft,  full  scale 

Separation  distances  of  the  store  eg  from  the  pylon-axis  system 
origin  in  the  positive  XP,  YP,  and  Zp  directions,  respectively,  ft, 
full  scale 

Translation  of  the-store  hook  from  the  carriage  position  in  the  Xp 
direction,  ft,  full  scale 

Physical  dimensions  of  the  CTS  support  sting,  in.  (see  Fig.  A-2) 


For  restricted  motion,  distance  store  must  travel  along  rail  in  a 
translate-only  mode,  ft,  full  scale 

For  restricted  motion,  distance  aft  hook  must  travel  along  rail 
before  becoming  free  of  rail,  ft,  full  scale 

Positions  of  the  CTS  pitch  center  with  respect  to  its  midpoint  of. 
travel,  positive  in  the  positive  Xt,  Yt,  and  Zt  directions, 
respectively,  in. 

Distances  (excluding  deflections)  from  the  CTS  pitch  center  to  the 
store  eg,  positive  in  the  positive  Xt,  Yt,  and  Zt  directions, 
respectively,  in. 

Positions  of  the  store  eg  with  respect  to  the  reference-axis  system 
origin  in  the  Xref,  Yref,  and  Zref  directions,  respectively,  ft,  full 
scale 

Separation  distances  of  the  store  eg  from  the  flight-axis  system 
origin,  positive  in  the  positive  Xt,  Yt,  and  Zt  directions, 
respectively,  ft,  full  scale 
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Xtp,o>Ytp,0 

Ztp.o 


X,,Y,.Z, 
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Y™,Z, 


cg>^cg 


AYC,AZC 

Ym  ,  Zm 


YTp 


Distances  from  the  midpoint  of  CTS  pitch  center  travel  to  the  store 
eg  at  carriage,  positive  in  the  positive  Xt,  Yt,  and  Z,  directions, 
respectively 

Distances  from  the  lanyard  attachment  point  on  the  aircraft  to  the 
store  eg  at  carriage,  positive  in  the  positive  XB,  YB,  and  ZB 
directions,  respectively,  ft,  full  scale 

Axial  distances  from  the  balance  electrical  center  to  the  model 
mass  center  in  the  pitch  and  yaw  planes,  respectively,  in.,  positive 
in  the  positive  XB  direction 

Full-scale  lateral  and  vertical  distances  from  the  balance  axis  to  the 
store  eg,  respectively,  ft,  positive  in  the  positive  YB  and  ZB 
directions 

Input  constants  used  to  modify  touch  point  coordinates  in  the 
lateral  and  vertical  directions,  respectively,  in. 

Lateral  and  vertical  distances  from  the  balance  axis  to  the  store  eg, 
respectively,  in.,  model  scale,  positive  in  the  positive  YB  and  ZB 
directions 

Touch  point  value  in  the  lateral  direction,  counts 


y,z  Lateral  and  vertical  distances  from  the  balance  electrical  center  to 

the  model  mass  center,  respectively,  in.,  positive  in  the  positive  YB 
and  ZB  directions 

ZERO  SET  Run/point  number  of  the  air-off  set  of  instrument  readings  used  in 

data  reduction 


ZelZk 


Zl 


Input  values  of  stroke  length  (or  time  of  action)  of  the  forward 
and  aft  ejectors,  respectively,  sec  or  ft 

Input  value  of  lanyard  length,  ft 


Zl,c 


ZSTEP 


Zt,c 


Current  value  of  straight-line  distance  between  lanyard  attachment 
points  on  store  and  aircraft,  ft 

Input  value  of  lanyard  length  at  which  the  data  cycle/integration 
time  interval  increase  option  is  exercised,  ft 

Thrust  module  event  designation  parameter,  ft  (see  Appendix  F) 
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Ejector  module  event  designation  parameters,  sec  or  ft  (see 
Appendix  G) 

Current  values  of  stroke  length  (or  time  of  action)  of  forward  and 
aft  ejectors,  respectively,  sec  or  ft 

Aircraft  model  angle  of  attack  and  sideslip  angle  relative  to  the 
free-stream  velocity  vector,  respectively,  deg 

Store  model  angle  of  attack  and  sideslip  angle  relative  to  the  free- 
stream  velocity  vector,  respectively,  deg 

Aircraft  angle-of-attack  value  recorded  in  the  touch  point  file,  deg 

Sidewash,  upwash,  and  crossflow  angles  with  respect  to  the  store 
longitudinal  axis,  respectively  positive  inboard  (left  wing),  up,  and 
clockwise  (from  -Zb  axis)  as  seen  by  the  pilot,  deg 

Simulated  aircraft  dive  angle,  positive  for  decreasing  altitude,  deg 

Indicated  angle  (in  pitch;  calculated  using  Cp>e)  between  the 
projection  of  the  local  flow  velocity  vector  onto  the  probe  Xb-Zb 
plane  and  the  probe  Xb  axis,  positive  for  a  velocity  vector 
component  in  the  negative  Zb  direction,  deg 

Angle  between  the  store  longitudinal  axis  and  its  projection  in  the 
Xp-Zp  plane,  positive  when  the  store  nose  is  to  the  right  as  seen  by 
the  pilot,  deg 

Angle  between  the  store  longitudinal  axis  and  its  projection  in  the 
Xp-Zp  plane,  positive  when  the  store  nose  is  to  the  right  as  seen  by 
the  pilot,  deg 

Angle  between  the  store  longitudinal  axis  and  its  projection  in  the 
XpZi  plane,  positive  when  the  store  nose  is  to  the  right  as  seen  by 
the  pilot,  deg 

Initial  input  angle  between  the  store  longitudinal  axis  and  its 
projection  in  the  Xj-Zj  plane,  positive  when  the  store  nose  is  to  the 
right  as  seen  by  the  pilot,  deg  (postlaunch  only) 

Initial  calculated  angle  between  the  store  longitudinal  axis  and  its 
projection  in  the  XpZi  plane,  positive  when  the  store  nose  is  to  the 
right  as  seen  by  the  pilot,  deg 
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?7r  CTS  rig  yaw  angle,  deg 

7]t  Angle  between  the  store  longitudinal  axis  and  its  projection  in  the 

Xt-Zt  plane,  positive  when  the  store  nose  is  to  the  right  as  seen  by 
the  pilot,  deg 

d  Angle  between  the  store  longitudinal  axis  and  its  projection  in  the 

XF-YF  plane,  positive  when  the  store  nose  is  raised  as  seen  by  the 
pilot,  deg 

Ad  Angle  between  the  store  longitudinal  axis  and  its  projection  in  the 

Xp-Yp  plane,  positive  when  the  store  nose  is  raised  as  seen  by  the 
pilot,  deg 

A0ck  Value  of  A 6  during  pivot  restrained  motion,  deg  (see  Appendix  I) 

6\  Angle  between  the  store  longitudinal  axis  and  its  projection  in  the 

X[-Y[  plane,  positive  when  store  nose  is  raised  as  seen  by  the  pilot, 
deg 

6P  Rotation  angle  of  the  aircraft  (flight-axis  system)  during  a 

pullup/pushover  maneuver,  deg  (see  Fig.  10) 

A0r  For  restricted  motion,  pitch  angle  through  which  store  must  rotate 

before  release,  deg 

d s  .  Prebend  angle  of  the  CTS  support  sting  in  the  pitch  plane,  deg 

dj  .  Angle  between  the  local  flow  velocity  vector  and  the  negative  Xp 

axis,  deg 

X  ■  Aircraft  model  scale  factor 

v  Angle  between  the  projection  of  the  store  longitudinal  axis  in  the 

XF-ZF  plane  and  the  XF  axis,  positive  when  the  store  nose  is  raised 
as  seen  by  the  pilot,  deg 

Av  Angle  between  the  projection  of  the  store  longitudinal  axis  in  the 

Xp-Zp  plane  and  the  Xp  axis,  positive  when  the  store  nose  is  raised 
as  seen  by  the  pilot,  deg 

v i  Angle  between  the  projection  of  the  store  longitudinal  axis  in  the 

Xi-Zi  plane  and  the  Xj  axis,  positive  when  the  store  nose  is  raised 
as  seen  by  the  pilot,  deg 
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Initial  input  angle  between  the  projection  of  the  store  longitudinal 
axis  in  the  Xj-Z]  plane  and  the  Xi  axis,  positive  when  the  store  nose 
is  raised  as  seen  by  the  pilot,  deg  (postlaunch  only) 

Initial  calculated  angle  between  the  projection  of  the  store 
longitudinal  axis  in  the  Xj-Zj  plane  and  the  Xj  axis,  positive  when 
the  store  nose  is  raised  as  seen  by  the  pilot,  deg 

CTS  rig  pitch  angle,  deg 

Angle,  between  the  projection  of  the  store  longitudinal  axis  in  the 
Xt-Zt  plane  and  the  Xt  axis,  positive  when  the  store  nose  is  raised 
as  seen  by  the  pilot,  deg 

Density  at  the  simulated  altitude,  slug/ft3 

Indicated  angle  (in  yaw;  calculated  using  CPi(J)  between  the 
projection  of  the  local  flow  velocity  vector  onto  the  probe  XB-YB 
plane  and  the  probe  XB  axis,  positive  for  a  velocity- vector 
component  in  the  positive  YB  direction,  deg 

Angle  between  the  store  lateral  (YB)  axis  and  the  intersection  of  the 
YB-ZB  and  Xp-Yp  planes,  positive  clockwise  looking  upstream,  deg 

Angle  between  the  store  lateral  (YB)  axis  and  the  intersection  of  the 
Yb-Zb  and  Xp-YP  planes,  positive  clockwise  looking  upstream, 
deg 

Simulated  aircraft  bank  (roll)  angle,  positive  clockwise  looking 
upstream,  deg 

Angle  between  the  store  lateral  (YB)  axis  and  the  intersection  of  the 
Yb-Zb  and  Xi-Yi  planes,  positive  clockwise  looking  upstream,  deg 

Separation  model  roll  orientation  with  respect  to  the  XB-ZB  plane, 
positive  clockwise  looking  upstream,  deg 

For  nonrolling  sting  applications,  the  calculated  angle  between  the 
true  and  simulated  roll  orientation  in  the  Euler  sequence,  deg  (see 
Appendix  L) 

Angle  between  the  projection  of  the  store  longitudinal  axis  in  the 
Xp-Yp  plane  and  the  Xp  axis,  positive  when  the  store  nose  is  to  the 
right  as  seen  by  the  pilot,  deg 
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A  \fr 

Angle  between  the  projection  of  the  store  longitudinal  axis  in  the 
Xp-Yp  plane  and  the  Xp  axis,  positive  for  store  nose  to  the  right  as 
seen  by  the  pilot,  deg 

h 

Angle  between  the  projection  of  the  store  longitudinal  axis  in  the 
Xi-Yi  plane  and  the  Xj  axis,  positive  for  store  nose  to  the  right  as 
seen  by  the  pilot,  deg 

’/'S 

Prebend  angle  of  the  CTS  support  sting  in  the  yaw  plane,  deg 

CO 

Angle  between  the  store  vertical  (ZB)  axis  and  the  intersection  of 
the  Yb-Zb  and  Xp-Zp  planes,  positive  for  clockwise  rotation  when 
looking  upstream,  deg 

Aco 

Angle  between  the  store  vertical  (Zb)  axis  and  the  intersection  of 
the  Yb-Zb  and  Xp-Zp  planes,  positive  for  clockwise  rotation  when 
looking  upstream,  deg 

Wl 

Angle  between  the  store  vertical  (ZB)  axis  and  the  intersection  of 
the  Yb-Zb  and  Xi-Z]  planes,  positive  clockwise  looking  upstream, 
deg 

WI,0 

Initial  input  angle  between  the  store  vertical  (ZB)  axis  and  the 
intersection  of  the  Yb-Zb  and  Xi-Zj  planes,  positive  for  clockwise 
rotation  when  looking  upstream,  deg  (postlaunch  only) 

Angle  of  the  simulated  ejector  force  with  respect  to  the  store  XB- 
,  ZB  plane,  deg 

C0o 

Initial  calculated  angle  between  the  store  vertical  (Zb)  axis  and  the 
intersection  of  the  Yb-Zb  and  Xi-Zi  planes,  positive  for  clockwise 
rotation  when  looking  upstream,  deg 

COR 

CTS  rig  roll  angle,  deg 

Acoxr 

For  nonrolling  sting  applications,  the  calculated  angle  between  the 
true  and  simulated  roll  orientation  for  a  pitch,  yaw,  roll  sequence, 
deg  (see  Appendix  L) 

cot 

Angle  between  the  store  vertical  (Zb)  axis  and  the  intersection  of 
the  Yb-Zb  and  XrZt  planes,  positive  for  clockwise  rotation  when 
looking  upstream,  deg 

(  )j 

Represents  the  ith  point 
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(  ) 

A  single  dot  denotes  the  first  derivative  of  a  parameter  with  respect 
to  time 

INERTIAL-AXIS  SYSTEM  DEFINITIONS 

Coordinate  Directions 

j 

x, 

Parallel  to  the  aircraft  flight  path  direction  at  store  release, 
positive  forward  as  seen  by  the  pilot 

Yi 

Perpendicular  to  the  Xj  and  Z\  directions,  positive  to  the  right  as 
seen  by  the  pilot 

z, 

Parallel  to  the. aircraft  plane  of  symmetry  and  perpendicular  to  the 
aircraft  flight  path  direction  at  store  release,  positive  downward  as 
seen  by  the  pilot 

Origin 

The  inertial-axis  system  origin  is  coincident  with  the  store  eg  at  release  and  translates 
along  the  initial  aircraft  flight  path  direction  at  the  free-stream  velocity.  The  coordinate  axes 
do  not  rotate  with  respect  to  the  initial  aircraft  flight  path  direction.  /  ' 

FLIGHT-AXIS  SYSTEM  DEFINITIONS 

Coordinate  Directions 


XF 

Parallel  to  the  current  aircraft  flight  path  direction,  positive 
forward  as  seen  by  the  pilot 

Yf 

Perpendicular  to  the  Xp  and  Zp  directions,  positive  to  the  right  as 
seen  by  the  pilot 

Zf 

Parallel  to  the  aircraft  plane  of  symmetry  and  perpendicular  to  the 
current  aircraft  flight  path  direction,  positive  downward  as  seen  by 
the  pilot 

Origin 

The  flight-axis  system  origin  is  coincident  with  the  store  eg  at  release.  The  origin  is  fixed 
with  respect  to  the  aircraft  and  thus  translates  along  the  current  aircraft  flight  path  at  the 
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free-stream  velocity.  The  coordinate  axes  rotate  to  maintain  alignment  of  the  XF  axis  with 
the  current  aircraft  flight  path  direction. 

PYLON-AXIS  SYSTEM  DEFINITIONS 

Coordinate  Directions 


Xp 

Parallel  to  the  store  longitudinal  axis  at  release  and  at  constant 
angular  orientation  with  respect  to  the  current  aircraft  flight  path 
direction,  positive  forward  as  seen  by  the  pilot 

Yp 

Perpendicular  to  the  Xp  direction  and  parallel  to  the  Xp-Yp  plane, 
positive  to  the  right  as  seen  by  the  pilot 

zP 

Perpendicular  to  the  Xp  and  Yp  directions,  positive  downward  as 
seen  by  the  pilot 

Origin 

The  pylon-axis  system  origin  is  coincident  with  the  flight-axis  system  origin  and  the  store 
eg  at  release.  It  is  fixed  with  respect  to  the  aircraft  and  thus  translates  along  the  current 
aircraft  flight  path  at  the  free-stream  velocity.  The  coordinate  axes  rotate  to  maintain 
constant  angular  orientation  with  respect  to  the  current  aircraft  flight  path  direction. 

STORE  BOD  Y-AXIS  SYSTEM  DEFINITIONS 

Coordinate  Directions 

XB  Parallel  to  the  store  longitudinal  axis,  positive  direction  is 

upstream  at  store  release 

Yb  Perpendicular  to  XB  and  ZB  directions,  positive  to  the  right 

looking  upstream  when  the  store  is  at  zero  yaw  and  roll  angles 

ZB  Perpendicular  to  the  XB  direction  and  parallel  to  the  aircraft  plane 

of  symmetry  when  the  store  and  aircraft  are  at  zero  yaw  and  roll 
angles,  positive  downward  as  seen  by  the  pilot  when  the  store  is  at 
zero  pitch  and  roll  angles 
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Origin 

The  store  body-axis  system  origin  is  coincident  with  the  store  eg  at  all  times.  The  XB,  YB) 
and  ZB  coordinate  axes  rotate  with  the  store  in  pitch,  yaw,  and  roll  so  that  mass  moments  of 
inertia  about  the  three  axes  are  not  time-varying  quantities. 

TUNNEL-AXIS  SYSTEM  DEFINITIONS 

Coordinate  Directions 

Xt 
Y, 

Zt 

Origin 

The  tunnel-axis  system  origin  is  located  on  the  tunnel  centerline  at  the  midpoint  of  axial 
travel  of  the  CTS  pitch  center  (Tunnel  Station  133.26) 

AIRCRAFT-AXIS  SYSTEM  DEFINITIONS 

Coordinate  Directions 


Parallel  to  the  tunnel  centerline,  positive  direction  is  upstream 

Perpendicular  to  the  Xt  direction  and  parallel  to  the  tunnel  top  and 
bottom  walls,  positive  to  the  left  when  looking  upstream. 

Perpendicular  to  the  Xt  and  Yt  directions,  positive  up 


XA 

Parallel  to  the  aircraft  longitudinal  axis  at  store  release  and  at 
constant  angular  orientation  with  respect  to  the  current  aircraft 
flight  path  direction,  positive  forward  as  seen  by  the  pilot 

Ya 

Perpendicular  to  the  XA  direction  and  parallel  to  the  Xp-Yp 
plane,  positive  to  the  right  as  seen  by  the  pilot 

Za 

Perpendicular  to  the  Xa  and  YA  directions,  positive  downward  as 
seen  by  the  pilot 

Origin 

The  aircraft-axis  system  origin  is  coincident  with  the  flight-axis  system  origin  and  the 
store  eg  at  release.  It  is  fixed  with  respect  to  the  aircraft  and  thus  translates  along  the  current 
aircraft  flight  path  at  the  free-stream  velocity.  The  coordinate  axes  rotate  to  maintain 
constant  angular  orientation  with  respect  to  the  current  aircraft  flight  path  direction. 
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REFERENCE-AXIS  SYSTEM  DEFINITIONS 

The  reference-axis  system  is  a  right-hand,  orthogonal  coordinate  system  whose 
coordinate  directions  (Xref.  YREf,  and  Zref)  and  origin  may  be  arbitrarily  selected  on  a 
test-by-test  basis.  The  most  common  alignment  of  the  coordinate  directions  is  parallel  to  the 
pylon-  or  aircraft-axis  system  coordinate  directions,  and  the  most  common  origin  locations 
are  coincident  with  the  store  eg  at  carriage  or  at  the  aircraft  fuselage  station,  waterline,  and 
buttock  line  zero  location. 
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